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Abstract 

The principle of minimal energy, which has been set up in the 
preceding papers for systems of non-identical particles (e. g. positro- 
nium), is now generalized to include also identical particles. Since 
the latter kind of particles feels also the exchange forces (besides the 
usual electromagnetic forces), one has to deal with non-zero exchange 
potentials which render the theory nonlinear, according to the non- 
Abelian character of Relativistic Schrodinger Theory (RST). However, 
the present extension of the variational principle refers only to the lin- 
earized version of RST in order to keep the calculations sufficiently 
simple. It is also demonstrated that in RST the Dirac particles can 
occur in fermionic and bosonic quantum states; and the mathematical 
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and physical consistency of the variational principle is validated for 
both types of states (concretely the fermionic hydrogen state 2p3 and 

the bosonic positronium state 2-'^Pi). 



PACS Numbers: 03.65.Pm - Relativistic Wave Equations; 
03.65.Ge - Solutions of Wave Equations: Bound States; 
03.65.Sq - Semiclassical Theories and Applications; 03.75.B 
- Matter Waves 
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I Introduction and Survey of Results 



The present paper is intended to promote the fluid-dynamic approach 
to the quantum phenomena. The conventional quantum theory is generally 
conceived to be a probabilistic framework for the description of microscopic 
matter (preferrably atoms and molecules). In order that this probabilistic 
description be most successful, it is also thought that it must be equipped 
with some counterintuitive elements, e.g. "... the so-called quantum pos- 
tulate, which attributes to any atomic process an essential discontinuity, or 
rather individuality, completely foreign to the classical theories and symbol- 
ized by Planck's quantum of action. This postulate implies a renunciation 
as regards the causal space-time co-ordination of atomic processes" [U |2]. 
It seems that this viewpoint of the early Bohr has been generally accepted 
in the meantime, because it seems to be safely supported by the acausal 
features of those observations which violate Bell's inequalities [3]. Indeed, 
Bohr even anticipated the occurence of such acausal effects by saying, "T/ie 
very nature of the quantum theory thus forces us to regard the space-time co- 
ordination and the claim of causality, the union of which characterizes the 
classical theories, as complementary but exclusive features of the description 
..." [1112]. Bohr's original idea of the impossibility of describing the quantum 
phenomena along the lines of a classical causal field theory (such as, e.g., 
classical electrodynamics) is held true up to the present days, as the following 
quotation may perhaps demonstrate most clearly: "// the classical concept 
of the space-time continuum were accepted, then quantum theory could not 
be considered complete, i. e., if it were accepted that the persisting objects of 
nature literally reside in a space-time continuum, ... then a complete scien- 
tific account of atomic phenomena would ... be required to describe whatever 
it was that is located at the points or infinitesimal regions of that continuum. 
Quantum theory does not do this, and hence a claim of completeness would 
be an abuse of language." [1] 

It seems to us that such assertions of Bohr, Stapp (and others) can dis- 
play only half of the truth. The point here is that, for the description of 
the quantum phenomena, we have to resort to a certain kind of dichotomic 
thinking which is usually termed '^wave-particle duality" [5J. By properly 
regarding this fact, one should feel being forced to admit that those claims 
about the unfeasibility of an ordinary space-time description for the quan- 
tum phenomena can refer only to the probabilistic (i. e. particle) approach, 
but by no means to the fluid-dynamic (i. e. wave) approach. Indeed, the 
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existence of a very successful fluid-dynamic approach (i. e. density functional 
theory [SIE!) demonstrates that the fluid-dynamic picture of quantum matter 
is at least as useful in atomic and molecular physics as the probabilistic view- 
point. Clearly, such an epistemological situation should bring forth sufficient 
motivation in order to further promote the fluid-dynamic approach which 
sometimes seems to have been overrun by the probabilistic way of thinking. 

In this sense, the subsequent discussion aims at the advancement of an 
alternative fluid-dynamic approach, i. e. Relativistic Schrodinger Theory 
(RST) [HI [9l [To]. In comparison to its competitors (i.e. the probabilistic 
approach and density functional theory), RST has obviously no problems at 
all with a completely relativistic formulation which, on the other hand, causes 
certain problems for both the conventional probabilistic approach and den- 
sity functional theory. Indeed it has already been demonstrated at various 
occasions that the Whitney sum construction (as opposed to the conventional 
tensor product construction) is perfectly compatible with the relativity prin- 
ciple (both special and general). This pleasant feature of RST is further 
worked out subsequently at many places, namely by passing over from any 
relativistic result to its non-relativistic approximation. Furthermore it will 
also be very instructive to explicitly verify the numerical coincidence of the 
RST and conventional predictions for some standard problem (i. e. the hydro- 
gen atom) where both the conventional theory and RST admit exact analytic 
solutions. The meaning of this is to demonstrate that RST owns at least the 
same degree of mathematical and physical consistency as the conventional 
quantum theory or density functional theory. 

However, the central point of the present investigation refers to the RST 
principle of minimal energy [9[ [TOl [11] . This is a variational principle on the 
basis of an appropriately defined energy functional i?[T] (both relativistic and 
non-relativistic) so that for the stationary bound systems the coupled set of 
mass eigenvalue and gauge field equations turns out as the extremal equa- 
tions due to that energy functional i?[T]. In some of the preceding papers, the 
principle of minimal energy has been set up for systems of non-identical par- 
ticles (e. g. positronium) where the exchange fields do trivially vanish. But 
of course one wishes to apply such a variational principle also to systems of 
identical particles which in RST are subjected to the action of the exchange 
forces. Since the latter are described in RST by additional gauge potentials, 
one would like to see the principle of minimal energy being validated also for 
this situation where the additional exchange potentials (B^) are non-trivially 
present besides the usual electromagnetic potentials (A^). The main result 
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of the paper is now that the inclusion of the exchange interactions is possible, 
but here the desired energy functional i?[T] must receive a slight modifica- 
tion and thus is no longer immediately equal to the spatial integral of the 
energy density Too (as the time component of the energy-momentum den- 
sity). Because of the considerable complications caused by the non-Abelian 
(and therefore nonlinear) character of RST, we present the proof only for the 
linearized version of the theory; but there should be no doubt that the non- 
Abelian generalization of the constructed energy functional i?[T] will exist 
also for the original non-Abelian theory. 

As a concrete demonstration of the usefulness and consistency of the 
proposed RST energy functional i?[T], we explicitly compute with its help 
the energy of both a fermionic and a bosonic field configuration (i. e. the 
hydrogen state 2p3 and the excited positronium state 2^Pi). 

These results are worked out in the following order: 



A. Linearization 

In sect.|lll the basic concepts of RST are briefly sketched in order to have 
the subsequent discussions sufficiently self-contained: the kinematics of the 
gauge fields is presented as the first part of the RST action principle whose 
second part consists then of the field equation for the pure states \1/ (for the 
mixtures see ref. [12] )• As in any field theory of the elemental particles, the 
field equations form the central part also for RST and they consist of the 
coupled system of matter equations and (non-Abelian) gauge field equations, 
see equations flll.31|l - flll.33|l below. Since these field equations emerge as the 
Euler-Lagrange equations due to the RST action principle, the corresponding 
Noether theorems automatically ensure the validity of the most important 
conservation laws, such as for current f lll.35[) and for energy-momentum 
Tfj^u ( III.45|l . The field-theoretic concept of energy- momentum is the starting 



point for the proper goal of this paper; namely to set up an energy functional 
(£^[T]) for the bound few-particle and many-particle systems in just such a 
way that the stationary form of the original Euler-Lagrange equations (due 
to the action principle) do appear simultaneously as the variational equations 
{''extremal equations") due to the desired energy functional E[t]^. 

Naturally, one thinks that the existence of an energy-momentum density 
Tfj^i, could provide a nearby basis for the definition of "energy" , i. e. one will 
first resort to the spatial integral of its time component 2oo(^)) see equa- 
tion flll.50[) below. Indeed, this procedure has been shown to be viable for 
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a system of non-identical particles (such as positronium [HI [IHl [H])- But 
for a system of identical particles the feasibihty of an analogous construction 
must be proven separately, because here the set of RST fields must be comple- 
mented by an additional field degree of freedom, i. e. the exchange potential 
Bfj_ which is responsible for the exchange interactions occurring exclusively 
among the identical particles. On the other hand, the inclusion of the ex- 
change interactions equips the theory with a truly non-Abelian character and 
this is a considerable complication. For this reason, we restrict ourselves in 
this paper to the linear approximation of the gauge field equations. Indeed 
this linearization process keeps the exchange fields present in the theory but 
on the other hand hnearizes the field equations for the gauge potentials, see 
the d'Alembert equations (lII.58a| -( ]II.58b|) below, so that the residual degree 
of complexity remains manageable. But clearly, there can be no doubt that 
an energy functional i?[T] of the desired kind would exist also for the original 
non-Abelian theory. The present treatment of the linearized theory, however, 
will show that for the inclusion of the exchange interactions the original idea 
of "energy" as the spatial integral flll.50[) of the time component TQ^ir) must 
be slightly modified. The nature of this modification is just one of the main 
results of the present paper. 



B. Electromagnetic and Exchange Fields 

Surely, the linearization of a non-Abelian gauge theory represents a con- 
siderable truncation of the original mathematical framework, and therefore it 
becomes necessary to reassure that the residual linear formalism is still log- 
ically consistent. Therefore in sect. lIlTl the linearized gauge field equations 
are analysed in detail. Since the goal of the paper aims at the stationary 
bound systems, this analysis is carried out in three-vector notation which 
presumes a certain space-time splitting so that the wave functions appear as 
the usual products of an exponential time factor and a time-independent wave 
function il)a{T)., see equation (]III.9[) below. Accordingly, the electromagnetic 
fields become completely time-independent, whereas the exchange fields 

still appear as the product of an exponential time factor and a time- 
independent spatial function, see equation flIII.6a|) below. This difference, 
however, does not disturb the time-independent curl and source equations 
which couple the three- vector fields to the charge and current densities as in 
ordinary Maxwellian electromagnetism, see equations flIII.13a|) - fllI1.14b|) be- 



low. It is merely the appearance of additional charge and current densities. 
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i.e. the exchange densities h^, h, which complete the ordinary densities of 
the Maxwelhan theory. But what is important here, this refers to the Pois- 
son and exchange identities which guarantee the numerical equality between 
the interaction energy M^^'^'^'^h^ of the gauge field modes with the matter 
modes on the one hand and of the gauge field energies E^^'^\ Eq'^'^ on the 
other hand; see equations (IIII.33|) - flIII.36|) below. This formal equivalence of 
the electromagnetic and exchange interaction energies once more points at 
the fact that in RST the exchange phenomena among identical particles are 
treated on the basis of real forces, quite analogously as their electromagnetic 
counterparts. 



C. Exchange Polarization 



The latter feature of RST is also responsible for a certain complication of 
the matter subsystem, in comparison to conventional quantum theory. The 
point here is that in RST all the ordinary densities (of charge, current, polar- 
ization, etc.) are accompanied by their exchange counterparts which generate 
the corresponding field strengths according to the gauge field equations, see 
below the equations (1III.3[) - (]III.4|I for the ordinary electromagnetic case and 
f lIII.10a[) - (]III.10b[) for the exchange case. This forces us to take into account 
also the exchange polarization effects which then entails a somewhat more 
intricate form of the mass eigenvalue equations. This complication becomes 
especially evident for the non-relativistiv form of the eigenvalue equations, cf. 
the simple Hamiltonian form fllll.74|l - fllll.75|) for the different (i. e. positively 
charged) particle vs. the extended form (1111.76a|l - pil.78b|l of both identical 
(negatively charged) particles. 



D. Mass and Energy Functionals 

Concerning now the search for an appropriate energy functional i?[T] for 
the stationary bound systems, one can deduce a kind of precursor directly 
from the (relativistic or non-relativistic) eigenvalue equations: this is the 
(relativistic) mass functional M[t]C^ (1111.43)1 . or the (non-relativistic) Pauli 
energy functional E[p] flIlI.129a)l - flIII.129b)l as the non-relativistic limit of 



M[T]C^. The point with these mass (or energy, resp.) functionals is that 
they do generate the corresponding eigenvalue equations (both relativistic 
and non-relativistic) via their variational (i.e. "extremal") equations. In- 
deed, the non-relativistic version may be viewed as the RST counterpart of 
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the well-known Ritz variational principle fHX [T3j if this is combined with 
the Hartree and Hartree-Fock approximations (see the discussion of this in 
ref. [H]). However, with respect to our original intention these relativistic 
mass and non- relativistic energy functionals M[t]C^ and Ep] suffer from a se- 
rious deficiency, namely it is not possible to deduce from them the gauge field 
equations! Therefore one has to undertake a fresh attempt from a somewhat 
different side. 



E. Exchange Interactions 

The new point of departure in sect. lIVl refers to the total energy i^T flll.50| 
as the integral of the total energy density '^'7oo(^) over all three-space. As it 
originally stands, this quanity cannot be adopted as the desired energy 
functional because the stationary field equations do not appear as the cor- 
responding variational equations. However, when the Poisson and exchange 
identities are added as constraints, one arrives at a modified form of E^, 
i.e. (1IV.12| . which almost meets with the expectations. Namely, the 

numerical value of that preliminary proposal E^^^ upon the solutions of the 
RST field equations is the same as for the original £^[t] (1II.50|) because both 
functionals E^J^-^ and i?[T] differ only by the Poisson and exchange identities 
and these are automatically satisfied by any solution of the field equations. 
But despite this pleasant feature of numerical equality of the latter pro- 



posal E^^i (]IV.12[) and the original E'fT] (1II.50|) . the functional i?™, cannot 



be accepted as the true solution of the problem if some of the particles are 
identical and thus undergo the exchange interactions. The reason for this 
rejection is two-fold: 

(i) The proposal -E^ij flIV.12|) contains not only the gauge field strengths 
but also the exchange potential B^. The usual expectation in field the- 
ory is that the energy concentrated in some field configuration should 
depend exclusively upon the field strengths (curvature components) 
and not upon the gauge potentials (connection components). 

(ii) The RST field equations do appear as the variational equations due to 
the proposed functional E^^^ (11V. 12|) only if no identical particles are 
present and thus the exchange potential is zero. Consequently, in 
the presence of identical particles {B^ ^ 0) the tentative proposal E^^!^-^ 
is not extremalized by the solutions of the RST eigenvalue problem! 
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In view of these problems with the right definition of the RST energy E'jx] 
we obviously have to abandon the energy concept (111.50| which is essentially 
based upon the existence of an energy-momentum tensor But fortu- 
nately, this is to be done in a minimal way: it is merely necessary to drop 
that questionable potential term from the proposal flIV.12|) so that the final 
result is represented by the energy functional E^,^1 (]IV.13|) . Indeed, this func- 
tional in both its relativistic form i?[x] (1IV-14|) and non-relativistic version 
E[T] (]IV.21|) has the corresponding RST field equations (matter plus gauge 
field) as its extremal equations and can therefore be viewed as the proper 
physical energy to be associated with any solution of the RST eigenvalue 
problem. But unfortunately, this reasonable definition of energy does then 
no longer coincide with the original RST energy E't (1II.50| which rests on 
the basis of the energy- momentum tensor 2^,^. Notice that the numerical 
equivalence of our final result E^:^'! and of the starting point i?T f!II.50|l is 
secured for all those situations where the exchange interactions are inactive 
non-identical particles). 



F. Application: Hydrogen State 2p3 

As a concrete application of the constructed energy functional, we con- 
sider in sect.|V]a one-particle system (the hydrogen state 2p3) and a two- 
particle system (the positronium state 2^Pi). Here it should not come as a 
surprise that for the one-particle states the RST predictions for the binding 
energy exactly coincide with their conventional counterparts. In fact, it has 
already been demonstrated for the one-particle systems by means of rather 
general arguments [151 IB] that the conventional one-particle energy (£fconv) 
always agrees with the corresponding total energy {E-r) of RST. Further- 
more, the latter always equals the mass eigenvalue (M2C^) so that for the 
one-particle systems one always has -Etonv = -^t = M2C^, see equation flV.53[) 
below. But despite this somewhat simple result it is nevertheless very instruc- 
tive to see in detail how the relatively complicated RST energy functional 
£'[T] manages to exactly coincide with the mass eigenvalue M2C^ on the one 
hand and with the conventional prediction -Etonv on the other hand. But 
beyond this numerical coincidence, the RST functional i?[T] provides a much 
more detailed picture of the intra-atomic situation. 
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G. Positronium State 2^Pi 



However, a more rigorous scrutiny of the RST energy functional i?[T] con- 
cerns the few-particle (or even many-particle) systems, because their treat- 
ment is on principle different in RST and the conventional theory (i. e. Whit- 
ney sum of one-particle bundles vs. tensor product of one-particle Hilbert 
spaces). As a test case we choose here positronium although its constituents 
(i. e. electron and positron) are non-identical particles and therefore the ex- 
change potential is zero. But nevertheless this is an absolutely non-trivial 
situation for testing the RST energy functional, since the bosonic character 
of positronium is in RST transferred to its constituents so that both the 
positron and the electron do occur in bosonic states (though being described 
by Dirac four-spinors) . Restricting ourselves here again to te non-relativistic 
limit, we can demonstrate that the RST prediction for the excited state 
approaches the conventional prediction (and the experimental value) up to 
a magnitude of (roughly) 10%, although we resorted to the spherically sym- 
metric approximation and to the simplest possible trial function with only 
two variational parameters, see equation flV.125[) below and the diagram on 
p. 11071 Since this magnitude of deviation is the same as for the analogous 
calculation of the groundstate (l"^So) energy in the preceding paper [10], we 
conclude that the use of the spherically symmetric approximation in principle 
generates such a 10-percent deviation which, however, would vanish through 
the use of better approximation techniques. 

Therefore it seems to us that the essential logical elements of RST (i. e. 
Whitney sum, energy functional and bosonic quantum states for Dirac par- 
ticles) must have something to do with physical reality. In this sense it 
would also be very interesting to develop more accurate (relativistic and 
non-relativistic) approximation techniques in order to check how close the 
RST predictions do really approach their conventional counterparts (and the 
corresponding experimental values). The outcome of such a competition be- 
tween a fluid-dynamic and a probabilistic view of the quantum phenomena 
would surely shed new light on the microscopic world. 

II Action Principle 

Once the basic field equations for RST had been set up, it was oberved 
that they may also be deduced from an action principle [12]. Indeed, the 
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existence of an action principle is always a very nice feature of a field theory 
because the relationship between the conservation laws and both the external 
and internal symmetries of the theory become elucidated in an elegant way 
(Noether theorem). But for the present case of RST, the consideration of 
the action principle is especially instructive also from another point of view, 
namely concerning its relationship with the principle of minimal energy. The 
latter principle is not only a pleasant feature of the general logical structure 
of the theory, but rather it is also of considerable practical relevance when 
one wants to develop variational techniques for the approximate calculation 
of the energy levels of bound systems [10] . 



1. Relativistic Field Kinematics 

The crucial element of any action principle, i. e. more concretely for RST 

SWnsT = 0, (II.l) 

is the Lagrangean density (£rst) which builds up the action integral M^st 
in the following way 

WnsT = Jd^xCnsT[^,A^]. (II.2) 

Since the RST field system is a coupled set of matter (\l/) and gauge fields 
(A^) , the corresponding Lagrangean appears naturally as the sum of a matter 
part (£d) and a gauge field part (Cq): 

CnsT[^,A^] = Cnm+CclA^]. (II.3) 

This splitting then induces a corresponding partitioning of the action integral 
Wrst (1II.2[) into two contributions, i. e. 

Wkst = Wd + Wg, (II.4) 

with the self-evident arrangement 

Wg[A^] 



(II.5a) 
(II.5b) 
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Here, the matter Lagrangean Cd looks as follows: 



tfic r — — 1 — 

c^[^] = — ^^r^(p^$) - (p^^)r'^J - ^Mc^"^ , (11.6) 

where the A^-particle wave function \I/ is the direct (i. e. Whitney) sum of 
the single-particle wave functions ijja (a = ^, ■■■ ,N) 

= tjji®ij2® --^^N ■ (11.7) 

Thus, from the mathematical point of view, the wave function \E'(a;) is a 
section of a complex vector bundle over space-time with typical fibre 
where the individual wave functions ipa are chosen as Dirac's four-spinors. 
Furthermore, the total velocity operator is the direct sum of the ordinary 
Dirac matrices 7^^, 

r'^ = (-7^) ©7'^©7^© ... , (11.8) 

where the explicitly displayed combination of Dirac matrices would then refer 
to a three-particle system (A^ = 3), with the first particle (a = 1) carrying a 
positive charge unit (-w positron, proton, etc.) and the second (a = 2) and 
third particle (a = 3) being negatively charged (^ electrons), see ref. [T7] . 
The gauge covariant derivative of the wave function \E' is as usual 

V^^ = 9^^ + ^^^, (11.9) 

with the gauge potential Afj, [bundle connection) taking its value in the 
Lie algebra u(A^) of the structure group U{N). This structure group be- 
comes reduced to its subgroup U{N,N') C U{N), which is of dimension 
A^'^ + (A^ - A^') = ^ + N'iN' - 1) if A^' particles are identical and the other 
(A^ — A^') particles are different. For instance, for the three-particle sys- 
tem of two identical particles (A^' = 2, electrons) and one different particle 
(A^ — A^' = 1, positron) mentioned below equation (111.8)1 the original nine- 
dimensional structure group f/(3) for three identical particles [TT] becomes 
reduced to the five-dimensional product group U{l)xU (2), where the Abelian 
factor U{1) refers to the positively charged particle; and the remaining four- 
dimensional factor group U{2) refers to the two (identical) electrons. The 
point with this bundle reduction is that the identical particles do not only 
undergo the Abelian electromagnetic interactions (described by [/(I)) but are 
additionally subjected to the non-Abelian exchange interactions (described 
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by the non-Abelian U{2), see ref. [T7] for a more detailed discussion). Ac- 
cordingly, the + N'{N' - 1) generators r„ (a = 1, ... , A + A^'(A^' - 1)) of 
the reduced Lie algebra u(A^, A^') are subdivided into the set of commuting 
electromagnetic generators Ta {a = 1, ... , A^) and the remaining N'{N' — 1) 
exchange generators Xk {k = 1, ... , A'(A' — 1)) so that the decomposition of 
the bundle comnection with respect to such a special Lie algebra basis 
appears in the following form: 

N N'{N'-1) 

A, = 5^A;r,+ Yl K^^- ("-10) 

a=l k=l 

In particular, for the three-particle system (1II.8| considered here one has 

3 

A = J^A^^T. + s^x-^x^ (11-11) 

a=l 

where the two exchange generators (fc = 1, 2) are taken as anti-Hermitean 
conjugates. Of course, the curvature J^i. of A^ (-^ field strength) 

T^u = V^A,-V,A^ + [A^,A,] (11.12) 

does then admit a similar decomposition 

3 
a=l 

where the F"'^^ are the (real- valued) electromagnetic field strengths and the 
(complex-valued) G^iy plays the role of the exchange field strength. The 
proper (continuous) gauge group for such an arrangement is now taken as 
the A^-fold Abelian product U{1) x U{1) x ... x U{1), which is the maximal 
Abelian subgroup of the original structure group U{N). The permutation of 
identical particles is an exact symmetry of the theory, with the corresponding 
discrete transformations to be considered as elements of the reduced structure 
group U{N,N') [n|. 

The (Whitney) sum structure of the A^-particle wave function \1/ flll.7[) 
lets now appear the matter Lagrangean Cd flll.6| as an (ordinary) sum of A 
single-particle contributions, i. e. for A^ = 3: 

C^m = Co(^i) + Cd(2) + C^is) , (11.14) 
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with the individual terms being given by 



^D(l) = 


ihc ~ 


2~ . 


Cb(2) = 


ihc ' 


'Y . 


^D(3) = 


ihc ' 





(11.15a) 
(IL15b) 
(II.15c) 



Here, the gauge-covariant derivatives {D^ipa) of the single-particle wave func- 
tions ifja {oL = 1,2,3) arise from the three-particle derivative T)^^ flll.9|) via 
the sum structure of ^ flll.7|) as 



(11.16) 



(II. 17a) 
(11.17b) 

(II. 17c) 



(D^^i) © (D^^2) © , 
with the individual derivatives looking as follows 

D.^l^i = d.iJi - I [A% + A\] 
D^^2 = d^ij2 - I \A\ + Ay - iB^^^ 
D^^i = d^i:3 - i + A\] ifj^ - iB^ilj^ • 
Consequently, the matter Lagrangean £d [^] PI.14)I is composed of four parts 

c^m = /:M[*] + /:kin[^] + 4i*] + ^£'\^] > (11.18) 

and any single one of these parts describes an important feature of the matter 
subsystem: 

(i) The rest-mass term Cm equips the RST field configurations with a rest- 
mass density which implies a rest-mass renormalization (see below) 



Cum 



-MpC^ ifjll/jl - MeC^ 



^2^2 + V^3^3 



(11.19) 



(ii) The kinetic term £kin takes into account the four- dimensional kinetic 
energy of the particles and differs in sign for differently charged parti- 
cles 

ihc 



+ 



ihc 



+ 



ihc 



(11.20) 
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(iii) The third term in (1II.18|) describes the electromagnetic interac- 
tions of the particles and is given by 

= hc[A\{kr+k.r)+A\{-kr+k,^)+A\{-ki^+kr 

(11.21) 

where the Dirac currents ka^ (a = 1, 2, 3) are defined as usual, i. e. 



k 



afj, 



(11.22) 



(iv) The last contribution C^^^ to the matter Lagrangean Cd flll.lSp is a pe- 
culiarity of RST because it describes the exchange interactions among 
the identical particles (i. e. electrons; a = 2,3) by means of the non- 
Abelian part of the gauge field theory: 



he 



(11.23) 



Observe here that the exchange generators x, x (!II.11[) do not com- 
mute, in contrast to the electromagnetic generators (a = 1,2,3), 
and thus they give rise to the non-Abelian character of the theory 
which is encoded in the exchange current h^ of the identical particles 
(a = 2,3): 



h,. 



'3 • 



(11.24) 



Once the general structure of the matter Lagrangean Cd is clarified, one 
may inspect the gauge field Lagrangean Cq in a similar way: 

he 

{a, (3 = 1, ... ,5) . 

Here, Kap is a covariantly constant fibre metric of the Lie algebra bundle and 
invariant under the reduced structure group [7(3, 2). Therefore, if the self- 
interactions of the particles are neglected (see ref.s [H [12]), this symmetric 
fibre metric {Kap = Kpa) adopts a very simple form, namely 

i^l2 = Kl3 = i^^23 = i^45 = -1 (11.26) 
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with the other (independent) elements being zero. Accordingly, the gauge 
field part /^Gf-^/i] of the RST Lagrangean £rst (!II.3|) appears as the sum of 
the electromagnetic and exchange contributions 



Cn + Cr 



(11.27) 



where the electromagnetic part (Cr) refers to the (real-valued) gauge field 
modes F% fULTSjl . i.e. 



he 



l_ll/ ' n.ij ' ii.ij 



(11.28) 



and similarly the exchange contribution [Cc] is related to the (complex- 
valued) exchange modes G^i, 



(11.29) 



For the subsequent variational procedure it is important to note that 
the electromagnetic field strengths -F"^^ and the exchange field strengths 
G^p (-^ curvature components) are linked to the connection components 
A"-^, through the following relations, to be deduced from the abstract 
equation (111.12[) : 



F] 



^.A'^ + t 



XA\^-i 



BuBf^, 



ByB^ 



G^. = \B,-V,B^ + i [A\ - A\] - I [Al - A\] B, . 



(11.30a) 
(11.30b) 

(11.30c) 
(11.30d) 



Naturally, the first field strength F^^^^ (]11.30a[) does not couple to the exchange 
potential B^ since the first electromagnetic mode {A^^, F^^^) is due to the 
first (differently charged) particle. On the other hand, the second and third 
modes F"^^^, F^^j, are linked also to the exchange potential B^ since the second 
(a = 2) and third (a = 3) particle do interact also via the exchange force. 
The exchange mode {B^,G^y) flll.30d|) does of course exclusively couple to 
both identical electrons (a = 2,3), but not to the positively charged particle 
(a = l)! 
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2. Relativistic Field Equations 

After the Lagrangean is fixed now in full detail, the variational proce- 
dure fill. If can be carried out in order to yield the complete set of coupled 
field equations. More concretely, the extremalization of the action integral 
WjisT f!II-2|) with respect to the wave function produces the A^-particle 
Dirac equation [T7] 



(11.31) 



or, resp., in component form for the presently considered three-particle sys- 
tem (A^ = 3) 



(II.32a) 
(II.32b) 
(II.32c) 



Furthermore, the extremalization of the action integral with respect to the 
bundle connection Afj, yields the non-Abelian Maxwell equations for the bun- 
dle curvature J-",, 



fJ,U 



-Airias Jv 

^2 



(11.33) 



(as = T-) , 

he 

i. e. in component form for the three-particle systems 









(II.34a) 






- B^G^, 


= Anasjl 


(II.34b) 




B'-G^, 


— B^G^y 


= Anasjl 


(II.34c) 






-Fy = Alias gu. 


(II.34d) 



Here, the emergence of the Maxwell currents j"^ {a = 1, ... ,5) needs an 
explanation because, up to now, the theory contains only the electromagnetic 
Dirac currents kafi flll.22[) and the exchange current /i^ flll.24|l . Namely, the 
point with the Maxwell four-currents j*^ is that they should obey a certain 
conservation law, i. e. in operator form 



v^^x = , 



(11.35) 
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in order that the Maxwell equations (1II.33[) be consistent. Decomposing here 
the current operator J^^ with respect to the chosen basis {r^, x, x} lets appear 
just the Maxwell currents j 



J, 



1. e. 



a=l 



SO that the conservation law (1II.35[) reads in component form 



(11.36) 



(11.37) 



Thus, one becomes forced to define the Maxwell currents j*^ in terms of the 
matter fields ipa in such a way that these source equations (111.37|) do emerge as 
an implication of the matter dynamics (111.31)1 itself! But this compatibility 
requirement for the matter and gauge field dynamics can be met just by 
means of the fibre metric K^is fill. 25)1 whose covariant constancy 

D^K^p = (11.38) 

helps solving the problem: First, define the RST currents ja^ through 



Jafi 



0) 



(11.39) 



which by a suitable choice of the generators turn out as nothing else than 
a combination of Dirac and exchange currents ka^, [8j 



k2fj. + k^n 
-ki^ + k2^ 

-h„ . 



(11.40a) 
(11.40b) 
(11.40c) 
(II.40d) 

(11.40e) 



Next, observe that these RST currents (111.39| do obey the following 
source equation: 



D^Ja^ = 



(11.41) 
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namely as an immediate consequence of the matter dynamics (1II.31[) . And 
finally, define the Maxwell currents (111.36)1 in terms of the RST currents 
jafi flll.39|) just with the help of the fibre metric Kaf3, namely through 



K 



al3 



(11.42) 



so that the covariant constancy ( 111.38)1 of the fibre metric ensures the com- 
patibility of the desired source equations fjll.37| for the Maxwell currents j*^ 
with those for the RST currents ja^ f)II.41)) . In this sense, the fibre metric 
Kap works also as a ^^compatibility tensor"; and the Maxwell currents can 
appear now in terms of the Dirac and exchange currents as 

(IL43a) 
(IL43b) 
(II.43c) 

(II.43d) 

3^,1 = -9t^ = ^ -^27/.^3 • (II.43e) 

But once a logical path from the Dirac and exchange currents to the 
Maxwell currents has been established, the source equations fjll.37)l for the 
latter can be transcribed to the coresponding source equations for the first 
kind of currents: 





= hf, = 




f = 


-h^i = 




f = 

J fJ. 




= -^37/.^/'3 




= = 


= ^■iltJi^2 


* 

= -Qi. = 


-h^ = 


= -^27/.^3 



V^ki^ = 



2fi 



B^h,., - B^'h,. 



= 
= 
= 



(II.44a) 
(II.44b) 

(II.44c) 

(II.44d) 



Observe here again that it is eclusively the electromagnetic current ki^ fjll.44a)) 
of the first (positively charged) particle which does obey a rigorous conser- 
vation law, whereas the other two (identical) particles fjll.44b)) - fjll.44c| must 
satisfy a more complicated source equation because of their exchange inter- 
actions. 

Besides charge conservation there is a further important conservation law 
which appears here as the source equation for energy-momentum: 



(11.45) 
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Since the RST field configurations represent a coupled system of matter and 
gauge fields, it is natural to assume that the total energy-momentum density 
^'^^T^u will be built up by a matter part '^'T^jy and a gauge field part 



(T 



(11.46) 



Indeed, the matter part has been identified as [T7] 
ihc 



- (P.^)r^^ + - {V^^)T,^ (11.47) 



and the gauge field part by 

he 



47ras 



(11.48) 



The (local) conservation law (1II.45|I comes now about through the mutual 
annihilation of the sources of both energy-momentum densities, i. e. 



- V ^ hcF^^^jj; . (11.49) 



Obviously the sources of the partial densities and are just the 

well-known Lorentz forces in non-Abelian form. 

It should now appear self-suggesting that the definition of the total energy 
(E't) of an RST field configuration is to be based upon the time component 
''^'Tqo of the energy-momentum density '■'^^T^^u, i- e. 



dV'^^Tn 



00 • 



(11.50) 



But since the total density '^-'3^1. is the sum of a matter part and a gauge 
field part, cf. (1II.46| . the total energy (111.50)1 naturally breaks up in an 
analogous way 



Et = Ed + Eq , 
with the self-evident definitions 



Ed= /dV(°)To 



00 



Er 



dV^^Too. 



(11.51) 

(II.52a) 
(II.52b) 
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However, such an arrangement obviously requires the sphtting up of space- 
time into space and time, which will readily be carried through for all the 
RST objects in the context of the bound systems. But the key question 
for the subsequent discussion will refer to the proposed definition flll.50|l for 
the energy concept in RST. Actually, we will readily show how this defini- 
tion flll.50| must be somewhat modified when the exchange interactions are 
taken into account. 



3. Linear Approximation 

The non-Abelian character of the A^-particle theory clearly implies a con- 
siderable complication when one wishes to look for concrete solutions of 
the field equations, i. e. the Dirac equation flll.31|l and the Maxwell equa- 



tion flll.33|) . On the other hand, it seems reasonable to assume that some 
aspects of the RST solutions can be discussed (at least in a qualitative way) 
by neglecting the nonlinearities which invade the theory through just its 
non-Abelian character. Therefore it is a nearby idea to consider first the 
linearized theory, namely by calculating certain energy levels in the linear 
approximation and comparing the corresponding results to the observational 
data. Such a linearization does not mean that one has to dispense with the 
exchange effects, because the exchange potential persists as part of the 
theory and it is only its field equation which becomes linearized. In this way 
one may deal with the exchange effects in a first approximation. 

The procedure of linearization refers mainly to the gauge field subsystem 
and consists in neglecting certain nonlinear terms emerging in that subsys- 
tem. In this sense, one neglects the commutators for the curvature defini- 
tion flII.12[) and puts in a simplifying way 

, (11.53) 
which truncates the component version flll.30a|) -f lll.30d[) to (a = 1,2, 3) 



n^ = %^^-V.A; (II.54a) 
G^, = V^B,-W,B^. (II.54b) 

Correspondingly, the component form (]II.34a|l - (lII.34d[) of the Maxwell equa- 
tion (]II.33|) becomes truncated to (a = 1, 2, 3) 

V^F^^, = A'Ka^jl (II.55a) 
W^G^^ = Anas K ■ (II.55b) 
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But this implies immediately that the non-Abelian imprints of the source 
system (lII.44a[) - (]II.44d[) become omitted and all the currents do obey local 
conservation laws in the strict sense: 



V^fc,^ = (II.56a) 
V^h^ = . (IL56b) 

Thus the linearized gauge field system adopts a shape very similar to ordinary 
electrodynamics; and if one takes over also the Lorentz gauge conditions 

VM';, = (II.57a) 
V^B^ = , (II.57b) 



then the ordinary Maxwell equations flll.55a[) - flll.55b[) are transcribed to the 



four-potentials in form of the well-known d'Alembert equations (a = 1,2,3) 

DA''^ = 47iasf, (II.58a) 
DBf, = Anas ^ ■ (II.58b) 

Moreover, the fact that the gauge field equations can be deduced from 
an action principle is not affected at all by the linearization procedure. In- 
deed, it is merely necessary to substitute the truncated curvature components 
(III.54a[) - flII.54b[) into the gauge field Lagrangean flll.25)) : and then the varia- 



tional procedure yields the linearized field equations flll.58a| - flll.58b|) in place 



of the original nonlinear ones fill. 34a| - fill. 34d|l . On the other hand, the mat- 
ter dynamics (]II.31|l . or its component form (]II.32a[) -( ]II.32c|) . resp., is left 



completely unaffected by the linearization process. 

Subsequently we will generalize the RST principle of minimal energy [9l 
[TT] in order to include also the exchange interactions, albeit only in the linear 
approximation; and thereby we will also have to modify somewhat the naive 
definition of total energy (1II.50| . 



Ill Stationary Bound Systems 

One of the most important applications of RST will surely be in the field 
of atomic and molecular physics, where the calculation of the relativistic 
energy levels must naturally receive the main interest. Therefore it becomes 
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necessary to subject the stationary RST states to a more thorough inspection. 
But this requires to break up the Lorentz covariant objects into their space 
and time components, which presumes the selection of some special Lorentz 
frame so that the manifest Lorentz covariance of the theory is lost. But 
the advantage of this reformulation of the relativistic theory in terms of 
three-tensors (of various rank) is that one can more conveniently obtain the 
stationary solutions of the RST field equations, both for the relativistic and 
non-relativistic situations. 



1. Gauge Field Subsystem 

The simplest space-time splitting refers to the four-potentials A"-^, which 
for the stationary states become time-independent and thus appear in the 
following form: 

{(%(f);-l,(r)} (IILl) 
(a = 1,2,3). 

A similar arrangement does apply also to the Maxwell four-currents 



so that the d'Alembert equations fill. 58a)) become split up into the Poisson 
equations, for both the scalar potentials *"Mo(r) 

A(%(f) = -47ras%(r) (IILS) 

and the three- vector potentials Aa{r) 

AAa{f) = -AnasUr) ■ (111.4) 

Recall here that the standard solutions of these equations are formally given 
by 

(%(f) = as /d¥' 3^ (ni.5a) 
Mr) = as fd'f' . (in.5b) 
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But somewhat more complicated is the space-time sphtting of the ex- 
change potentials because they are left time-dependent, even for the sta- 
tionary states: 



exp 



ct 



am 



where the ''exchange radius" au is defined as 

ttM 



(M2 - Ms)c 



(IIL6a) 
(III.6b) 

(III.7) 



The reason for such a time-dependence is that the exchange potential B^{x) 
must be a solution of the d'Alembert equation fill. 58b)) with the space-time 
splitting of the exchange current h^{x) ( 1II.24|) appearing as 



h^{x) = exp 



ct 



Am -I 



■ Kir) 



(111.8a) 
(111.8b) 



This specific time- dependence of the exchange current is directly induced by 
the single-particle wave functions tfjaix), which adopt their stationary forms 
in terms of the mass eigenvalues Ma as 



exp 



(III.9) 



Accordingly, the d'Alembert equation (1II.58b[) for the exchange potential 
becomes not reduced to the ordinary Poisson equation as for the electromag- 
netic potentials A"^, cf. pil.3|l -f iril.4|l . but rather does appear in the following 
form: 



A5o(r) + ^ So(f) = -47ras K{r) 



AB(f] 



B{r) = -Anas h{r) 



(III. 10a) 
(Ill.lOb) 
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The formal solution of this type of equation is easily found to be the following: 

r ^(fO-COsfJ^^') 

5o(f) = as /dV „^ >„;^ ^ (Ill.lla) 



h{r') • cos 



B{r) = as /dV — — ^ , (Ill.llb) 



which for infinite exchange radius (cm — ^ oo, i. e. M2 M3) tends to the 
massless form of the electromagnetic potentials *"Mo(r), Aa{r} fllll.5a)l - fllll.5b l. 

It is true, the particle interactions are organized via the (electromagnetic 
and exchange) potentials which, according to the principle of minimal cou- 
pling, are entering the covariant derivatives -D^V'a of the wave functions ipa as 
shown by equations (111.17a|l - (111.17c[) . But nevertheless it is very instructive 
to glimpse also at the field strengths F"-^^, G^y. Their space-time splitting is 
given by 

(III. 12a) 
(III.12b) 
(III.12c) 
(III.12d) 

and thus the linearized Maxwell equations (]II.55a|) - (]II.55b[) do split up in 
three- vector form (a = 1, 2, 3) to the scalar equations for the electric fields 

V-^„ = 47ras% (111.13a) 
V-X = 47ras4, (111.13b) 

and to the curl equations for the magnetic fields 

V X = Anas] a (111.14a) 

W xY + —X = Anash . (III.14b) 
au 

There is a pleasant consistency check for these linearized (but still relativis- 
tic) field equations in three-tensor form; namely one may first link the field 



Ea = 






Ha = 






X = 


{x^} = 


{^Oj} 


Y = 


{Y^} = 
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strengths to the potentials in three-vector notation (cf. (111.54a)) - (111.54b|l for 
the corresponding relativistic hnk): 



Xir] 

Ha{f] 

Y{f] 



-V<%(f) 

-VBn(f) + —B(f] 



V X Aa{r) 

V X B(f) , 



(111.15a) 
(111.15b) 

(III.15c) 
(III.15d) 



and then one substitutes these three-vector field strengths into their source 
and curl equations (]III.13a)) - (]III.14b|) . In this way one actually recovers the 
Poisson equations (1III.3[) - (11I1.4)) for the electromagnetic potentials '"Mq, A^, 
but the exchange potentials Bq, B require an extra discussion. 

The point here is that the Lorentz gauge condition flll.57b|) must in three- 
tensor form be consistent with the linearized source and curl equations for 
those three-vector field strengths. In order to become convinced of this, 
substitute first the magnetic exchange field strength Y (]III.15d)) into the curl 
equation ( ]III.14b[) and find 



V(V • S(f)) - AB{f) V5o(f 



B(f] 



47ras/i(r)- (111.16) 



Obviously it becomes necessary here to select some gauge condition for the 
magnetic exchange potential B in order that the latter equation (1111.16)1 can 
coincide with the already fixed equation flIII.10b|) for B. This requirement 
implies that one has to put 



(V'S(f)) - 

whose simplest solution is of course 



Bo{r) 







au 



V'B-—Bo{f) = 







(III.17) 



(III.18) 



But this is nothing else than the Lorentz gauge condition ( ]II.57b)l for the 
four-potential 5^ when its general time-dependence (IIII.6a|) is respected. 

Next, contract the curl equation flIII.14b|l by the gradient operator (V), 
observe also the generally valid identity 



V-(Vxy(f)) = 



(III.19) 
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and thus arrive at 



— V • X(f) = V • h{f) , (III.20) 



which by means of the source equation (IIII.13b|l leads to 



— ho{f) = V'h{r). (III.21) 

But this is again nothing else than the conservation law flll.56b|l for the 
four-current when its general time- dependence flIII.8a| for the stationary 
states is respected. In this way it is evident that the gauge part of the RST 
field equations remains a completely consistent framework for the stationary 
systems after the process of linearization has been carried through. 

Moreover, the exchange subsystem of the gauge part admits an important 
conclusion for the matter subsystem. Namely, observing the general defini- 
tion flll.24|l for the exchange current /i^, one finds for its time component 



/io(r) 

ho{r) = V'2(r)7o^3(r) = ^|(r)^3(r) ■ (111.22) 

Therefore integrating equation (]III.21|) over all three-space, with observation 
of the exchange radius om ( IIII.7|) . yields the following important result: 

{M2- M3) jd^fijl{f)tp3{f) = -i-jh{r)'AS, (III.23) 

iSoo) 

where 5*00 is some infinitely distant two-surface. Assuming now that the 
matter fields i^aif) (a = 2,3) decay rapidly enough (e.g. exponentially) at 
spatial infinity (r — > 00), the exchange flux through this infinitely distant 
surface will vanish, 

j h{f) 'dS = , (III.24) 

(5oo) 

and thus we are left with the orthogonality condition 

'dVV'|(f)^/'3(f) = . (III.25) 



28 



But observe here that this conclusion can be drawn only for those states 
V'alO (fl = 2,3) which have different mass eigenvalues (M2 7^ M3). 

Further information about the gauge field subsystem can be obtained by 

multiplying through the equations (]III.10a[) - (]III.10b[) by BQ{f) and B{f)^ 

resp., with subsequent integration over all three-space. In this way one finds 
from the first equation (lIII.10a[) 



d f 



1 * 



— B^{r)B^{r) - {VBo{r)) • (VSo(r)) 



47ras ld^rho{r)Bo{f) . 

(III.26) 



This result says that the right-hand side, i. e. the exchange mass equivalent 
of electric type (M^^c^), must be a real number: 



he I d^f Bo(f) ■ ho(f) = he I d^r Bo(f) ■ hoir] 



Boif) ■ hoif) + Boif) ■ hoir) 



[m.27) 



A similar conclusion can be drawn also from the second equation (lIII.10b|) 

for the exchange mass equivalent (M®c^) of the magnetic type: 



he jd^r B{r)'h{r) = he j d^f B{f) • h{f] 
he 



B{f) • h{r) + B{r) • h{r] 



^111.28) 



which emerges in the following form: 
1 



d?rB{r) • B{f) - /d¥ (V x 5(f)) • (V x B{f)) 



(III.29) 



-47ras / d^f h{f) • B{f) . 



On the other hand, the exchange system deals also with purely imaginary 
integrals, e. g. 



dV5(f)-X(f) = - /dV5(f)-X(f) . 



(III.30) 
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In order to verify this one simply substitutes here the exchange field strength 
of the electric type (i.e. X{f)) from equation (lIILlSbf and obtains through 
integration by parts 



dVS(f)-X(f) 



Am 



dV 



B{f)'B{f)-Bo{f)-Bo{f] 



(III.31) 



/dVS^S, 

Am 



Obviously, the integral on the left-hand side vanishes either for om oo, i. e. 
when both electrons occupy states with identical mass eigenvalues (M2 = M3), 
or it vanishes also for the situation where the exchange interactions of the 
electric (-w Bq) and magnetic type (-^ B) are of equal magnitude 



d-^fB^B,, 



dV 



Bo{f)-Bn(f)-B(f)'B(f] 



. (III.32) 



However, the magnetic-type interactions are usually much smaller than their 
electric counterparts, which then implies that the requirement ( 1111.32)1 can in 
general not be satisfied. One will subsequently have to return to this point 
of omission of the integral (1III.32|) in connection with the RST principle of 
minimal energy. 

In this context, the exchange identities of the electric (1III.26|I and mag- 
netic type (1111.29)1 will subsequently play the role of constraints and are for 
this purpose recast into a more suggestive form, namely for the electric type 



(h) 



d f 



he 

47rQ!p 



(Vi?o(r))-(Vi?o(r)) - — i?o(r) ■ Bo{f 

Bo{f) ■ ho{r) + Bo{f) ■ ho{f 



(III.33) 



he 
Y 



and similarly for the magnetic type 



(g) 



dV 



he 

47rQ!s 
1 



(V X 5(f))«(V X B{f)) + (V-S(f)) • (V-S(r)) 



a 



B{ryB{r) 



he 

T 



B{fyh{f) + B{f)'h{f) 



. 



(III.34) 
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It should be evident that similar identities can also be deduced from the elec- 
tromagnetic Poisson equations (1III.3| and (1III.4[) . namely the electric Poisson 
identity P HI] 



he 



/ he 



(V(%).(V(%) + (V(%).(V%o) + (V%o)-(V(%)" 
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(1). 



= 
(III.35) 



and its magnetic counterpart 



(m) 



he 
47ras 

he 

~ ~2 



(Vxli).(Vxl2) + (Vx^).(Vx/3) + (Vx/3).(Vx/i) 

Alin + js) + ^2-(j3 + Ji) + Min + /a)] I = . 

(III.36) 



Finally concerning the strength of the exchange interactions, the equa- 
tions ( lIII.lla)) -( ]III.llb]l tell us that for vanishing overlap of both wave func- 
tions -02 and the exchange density /iq and exchange current h do also 
vanish and hence the exchange potentials So, B become zero. Therefore the 
exchange interactions can exist only for overlapping wave functions! But 
there is another influence on their strength, and this refers to the magnitude 
of the exchange radius om f lIII.7)) , which roughly measures the energy differ- 
ence (M2C^ — M^e^) of both particles. In order to get a rough impression of 
this effect, assume the exchange density /io(0 to be of the exponential form 



ho{r) 



Vo exp 



r 

ax 



(III.37) 



where the length parameter ax measures the spatial extension of the over- 
lap region. This assumption then yields for the electric exchange potential 
Bo (Ull.llall at the origin (r = 0) 



Bn 



r=0 



47rQ;s Vo 



Am + 



I? 



(III.38) 



Consequently, for small enough an energy difference (M2 ~ M3, om ^ «x) 
the exchange potential is strongest; and conversely, for increasing energy 
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difference (om — ^ 0) the exchange interaction tends to zero. Clearly, when 
the wave functions are well separated (i. e. ax 0) then the exchange effects 
do also vanish. Summarizing, the exchange effects become negligible for well 
separated particles of sufficiently different energies. 



2. Matter Subsystem 

Similar to the gauge field subsystem, the matter subsystem exhibits also 
some peculiar features which are worth to be inspected in greater detail. 
When the stationary form (]III.9[) of the wave functions ipa{x) is inserted 



in the general matter dynamics (111.32a[) - (111.32c| . one is led to the mass eigen- 
value equations [17J for the spatial parts 4'a{'^) of the three wave functions. 
Since the first particle (a = 1) carries a positive charge, it is different from 
the other two particles (i. e. the electrons) and therefore does not feel the 
exchange force. Consequently, its mass eigenvalue equation does not contain 
the exchange potentials BqIt) and B{f) but exclusively the electromagnetic 
potentials ^'Mo(r) and Aa{f) (a = 2,3) due to the other two (identical) par- 
ticles. Furthermore, the preceding splitting of the gauge objects into their 
space and time components suggests a similar splitting of the spatial part 
ipa{^) of the wave functions which then spoils the manifest Lorentz covari- 
ance but nevertheless preserves the relativistic character of the theory. Such 
a splitting of the wave functions refers to the fact that the Dirac four-spinors 
i/jair) may be conceived as the direct sum of Pauli two-spinors ''^'(/?±(r): 

^a(f) = <V(0©*V(0 (III.39) 
(a= 1,2,3) , 



so that the matter dynamics (]II.32a[) - (]II.32c|) transcribes to these Pauli spinors 
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as follows [T7] : 



(1) 



^(r) + ^3(r) •a%(r) 



±M„ - Ml 



- li(f) + y43(r 



•a%(f) -S(f)-a®^±(f^ 



M2 ± Me 



M3±Me 



c-%(f) 
(IIL40a) 



(III.40b) 



(IIL40c) 



Observe here again the crucial point with this three-particle system; 
namely that both eigenvalue equations flIII.40b[) and flIII.40c|l for the two 
identical particles are directly coupled via the exchange potentials Boij) and 
B{r), so that the second particle's equation ( ]III.40b)l contains also the third 
Pauli spinor ®v^±(r) and vice versa; whereas the first particle (lIII.40a|l couples 
only indirectly to the other two (identical) particles, namely via the gauge 
fields '"Mo, Aa (a = 2, 3) being generated by just those identical particles. 
Therefore it is natural to suppose that such a direct (i. e. exchange) coupling 
will cause a much more intimate connection between the identical particles 
than is the case with their indirect coupling to the first (different) particle. 
Clearly this supposition requires a more thorough discussion below. 

But the coupled system of mass eigenvalue equations (lIII.40a|l - (lIII.40c|) 
yields now a first hint on how to set up an energy functional so that these mass 
eigenvalue equations can appear as the corresponding variational equations. 
It is true, the sum (M[t]) of the mass eigenvalues Ma (to be understood as 
functionals M^a] of the Pauli spinors) 



Ml 



[T]( 



-M[i]c2 + M[2]c2 + M[3]c2 



(III.41) 



can not yet be taken as the desired energy functional ^'[t]; but nevertheless 
the mass eigenvalue equations actually are the variational equations of this 
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total mass functional M[t]C^ (1III.41|) . under the constraint of wave function 
normalization (a = 1, 2, 3) 



Nr 



'D{a) 



dV(%(f) - 1 



. (III.42) 



In order to recognize this prehminary success more clearly, resolve the 
mass eigenvalue equations for the mass eigenvalues Ma (or M^a], resp.) by 
contracting those equations with the appropriate Pauli spinors and subse- 
quent integration. Then add up the three mass functionals M[a]C^ and thus 
find the total mass functional as [9] 



M[T]C2 = ) ■ MpC^ + Zf2) ■ MeC^ + Z^^) ■ M^C^ + 2Tkin + 2 



- 2 



(III.43) 



Here, the mass renormalization factors Zu) are given by 



d'f 



= d^f'il)a{r)"ifJa{r) , 



^111.44) 



and the kinetic energy Ikin is the sum of the three single-particle contributions 

^in = ?kin{l) + ^in(2) + ^in(3) , (III. 45) 

i. e. for the first particle 
he 



2 



%l (f )a • V^V (f ) + %l {r)a • (^^ 



= /d'^r '?/^i(r)7 • V'?/'i(r) 



(III.46) 



and analogously for the two identical particles (a = 2, 3) 



kin(a) 



-2— /d^f 
2 



Hot 



Vl(r)^* V'"W(0 + *"V-|(0^* V'V(^) 



^111.47) 
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Furthermore, the exchange mass equivalents M^^c^ and M®c^ of the elec- 
tric and magnetic type have already been defined through equations (1III.27|) 
and fllll.28|) : and their electromagnetic counterparts M^'^^c^ and M^^^c^ are 
given by 



M(V = y /dV 



%(r) 



(%(f) + (%(f)l - (%(f) ■ [(%(f) + %(f 
%o(r) + ^%if 



(III.48a) 



he 

y 



dV |fci(f) • 



A2{f)+A3{f] 

Ai(f) + A2(f) 



Ai(f)+A3(f] 



(IIL48b) 



Now with these arrangements it is an easy exercise to verify that the mass 
eigenvalue equations flIII.40a|l - flIII.40c[) indeed do represent the variational 
equations due to the variational principle 



6{M[T]C^) = 



(III.49) 



Here, the functional M[t]C^ is nothing else than the original M[t]C^ (1III.43 I 
merely to be complemented by the normalization constraints (]III.42)) : 



M[T]C2 = M[T]c2 + 5^AD(a)Ai,(a) , 



(111.50) 



a=l 



where the AD(a) are the usual Lagrangean multipliers. The comparison of the 
variational equations (1III.49|I with the mass eigenvalue equations (]111.40a|l - 
(]111.40c|l yields then immediately the equality of the Lagrangean multipliers 
and the mass eigenvalues, i. e. 



Ad(i) = Mic" 

Ad(2) = -M2C^ 
Vd{3) 



(111.51a) 
(III.51b) 
(III.51c) 



But observe here that, despite of this partial success, one cannot be sat- 
isfied with adopting the mass functional M[t]C^ flIII.50[) as the desired energy 
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functional because this mass functional does not admit to deduce the gauge 
field equations (1III.3| - (1III.4[) and (lIII.10a| - (]III.10b|) from the variational prin- 
ciple flIII.49|l . The reason is that the functional M[t]C^ does not contain any 
derivatives of the gauge fields '°Mo, Aa, Bq, B; but this is indispensable for 
the purpose of deducing those Poisson equations for the gauge fields. Nev- 
ertheless, a certain non-relativistic limit of the variational approach flIII.49[) 
is sometimes used in conventional quantum theory (i. e. the Hartree and 
Hartree-Fock approximation), but this does not in every case yield acceptable 
results |Tlj. Thus the non-relativistic limit of RST requires a more thorough 
inspection, especially concerning the properties of the matter subsystem. 



3. Polarization of Matter 

Since the total mass functional M[t]C^ (]III.50[) does correctly reproduce 
the matter part of the RST dynamics, it can be expected to lend itself to an 
elucidation of how the matter interacts with the gauge fields. First, consider 
here the positively charged particle (a = 1), which formally contributes to the 
total mass functional M[t]C^ (IIII.43[) exclusively through its electromagnetic 

interactions M^^'^^^c^ 



Mj^^c'' 



y /dV-(% 

he 
~2 



— Id^fk 



y Id'rh 



A2 + A3 



(III.52) 



cf. equations flIII.48a)) - flIII.48b)l . Now observe here the Gordon decomposi- 
tion p2] of the first Dirac current kifj, into its convective and polarization 
parts, i. e. 

ht, = -qif, - • (III. 53) 

Here, the convection current ((/i^) is given in terms of the covariant deriva- 
tive (111.17a[) as 

ih 



9iM 



2MpC 



(III.54) 



while the polarization current (si^) is found to be of a divergence form 



(III.55) 
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with the electromagnetic polarization tensor S 
tide being defined through 



(1) 



S^u^n) of the first par- 



c(i) 



(III.56) 



Now by means of this decomposition flIII.53)) of the first Dirac current 
/ci^, the electromagnetic mass equivalent flIII.52|) can also be split up into 
two parts: 



he 



_|_ j^ti 



he 
Y 



d^fsi. 



^^2^ _|_ j^ilJ. 



^111.57) 



Obviously, the convective part Mj^jjc^ represents here the interaction of the 
purely translational degrees of freedom of the particle (-w point particle), 
whereas the polarization part Mj^c^ refers to the particle's dipole properties 
(of both the electric and magnetic type). The latter property is elucidated 
more clearly by observing the divergence form flIII.55| of the polarization 
current. This latter form admits namely to recast the polarization part m'^^c? 
of the mass equivalent (1III.57|I into the following shape: 



-— /d^fs 



J^IJ- _|_ 

E2 + 



(m) 



^^1 



#2 + i?3]} , (111.58) 



with the electric (e) and magnetic (m) polarization densities ^"^Si = {'"^S*/} 
and '™'S'i = {'""IS/} being related to the space and time components of the 
original four-tensor S^]]u ( 1III.56|) through 



S 



(1) 



ih 



2MpC 



2MpC 



(III.59a) 
(III.59b) 
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The time component *^%o of the polarization current Si^ (1III.55| is then noth- 
ing else than the source of the electric dipole density ^°^Si (1111. 59a| 



(1), 



)so = -V • % , (111.60) 



whereas the curl of the magnetic dipole density '"'iSi ( 1111.59b)l is related to 
the space part si = {s/} by 

si = V X (-"'S^i . (111.61) 

Thus the conclusion is that a Dirac particle owns both a translational 
(~ q^) and a rotational (~ s^) handle to be acted upon by the electromag- 
netic fields. This will become even clearer below in the course of looking for 
the non-relativistic approximations. 

However, the point with this polarization phenomenon is now that in 
RST the identical Dirac particles are subjected additionally to the exchange 
forces; and since the latter are described by the same gauge formalism as 
their electromagnetic counterparts, it is logical to assume that those exchange 
forces do also act upon both the translational and the rotational degrees of 
freedom of the identical particles. Indeed, one is easily convinced of this by 
inspecting now in more detail these exchange interactions of the second and 
third particle (a = 2, 3). Quite similarly to the electromagnetic case (1111.52)1 . 
the exchange contribution to the total mass functional Mjtjc^ ( 1111. 43[) can 
also be rewritten in terms of a Lorentz invariant, i.e. 



(111.62) 



cf. equations (]111.27jl -( !ITL28|l . 

Next, one works out the Gordon decomposition of the exchange current 
( lll.24[) in a way quite analogous to that for the electromagnetic current 



ki^i (1111. 53|l . namely by substituting the wave functions ■02 and il)^ from 



the second and third Dirac equation ( ]11.32b[) -( ]11.32c)l . This then yields in 
a straightforward manner the following form for the exchange current h^: 

K = K + [^'^ - ^'1 Z,. - iB^ [S^% - . (111.63) 

Of course, the first two terms are here just the exchange counterparts of the 
convection and spin currents qi^ and Si^ of the first particle ( 1111. 53)l . i.e. 
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more concretely 



ih 



2MeC 



fiu -I 



with the exchange polarization tensor Z^^, being given by 



ih 

'mTc 



(III.64a) 
(III.64b) 

(111.65) 



When this decomposition flIII.63|) of the exchange current /i^ is now 
substituted into the mass equivalent ( ]1I1.62|) . the latter splits up into four 
contributions, where the first two correspond to the translational and rota- 
tional degrees of freedom of the two identical particles; for the rotational 
case see the electromagnetic counterpart (]III.58|l where the electric and mag- 
netic fields Ea, Ha (a = 2, 3) have to be replaced by their exchange coun- 
terparts X, Y ( ]III.15b)l and ( ]III.15d)l . and similarly for the dipole densities. 
However, in addition to these exchange counterparts of the electromagnetic 
forces, there arise two new terms in the Gordon decomposition of the ex- 
change current h^ (1III.63[) . i.e. the terms with the polarization tensors Z^j^^y, 
S^%, and S^%. Therefore, if the exchange current h^ (1111. 63|) is inserted into 
the exchange mass equivalent ( 1III.62)I . there will emerge terms quadratic in 
the gauge potentials which then implies certain non-linearities in the the- 
ory (as it is to be expected for a non-Abelian gauge theory). Naturally, 
these nonlinearities will entail considerable computational complications so 
that one will omit them for the present linear approximation. Thus for our 
present linearized approach one can resort to the following truncated form of 
the exchange current h^ ( 1III.63 I 



h„ 



bf^ + 



(III.66) 



which appears then as the formal counterpart of the electromagnetic currents 
such as that for the first particle (1III.53|) . This will readily be elaborated in 
greater detail in connection with the non-relativistic limit. 



4. Non-Relativistic Limit 

The value of any physical theory must surely be measured by the ac- 
curacy of its predictions with respect to the observational data. But here 
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it is important to relate this question to the various degrees of accuracy 
which one wishes to demand from a theoretical formalism. Indeed, one of 
the most frequent gradations of accuracy refers to the relativistic vs. non- 
relativistic approaches. Naturally, one will prefer to inspect and test first the 
non-relativistic limit of the theory in question, before its fully relativistic con- 
sequences are considered. In the present context of the principle of minimal 
energy, this viewpoint would entail that one first has to deduce the non- 
relativistic limit of this principle from its original relativistic form, namely in 
order to test the corresponding non-relativistic predictions, e. g. in the field 
of atomic physics. But for this purpose, one will find it logically desirable 
that the non-relativistic form of the original RST field equations should be 
identical to those variational equations which may be deduced from the non- 
relativistic limit of the principle of minimal energy. Or in other words, the 
following commutative arrangement should be true: 



Relativistic 


non-relativistic ^ 


Non- Relativistic 


Field Equations 


approximation 


Field Equations 



Relativistic Principle 
of Minimal Energy 


non-relativistic ^ 


Non-Relativistic Principle 
of Minimal Energy 


approximation 



For the verification of this demand one may first turn to the mass eigen- 
value equations whose non-relativistic approximation is easily found by sim- 
ply eliminating the "negative" Pauli components (f) and thus concen- 
trating on the remaining eigenvalue equations for the "positive" components 
^"•^Lp^{f). For this purpose, one formally resolves the original mass eigenvalue 
equations (]III.40a|) - (]III.40c[) for the negative Pauli components (r) in the 
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following approximate way: 



(1) 



a ' 



(2), 



2MpC 

ih 



(1), 



(3), 



ih 



2MpC 



a ' 



V + iAi 

V + iAii 

V + iAiii 



V+(r) 



(2), 



h 



2MeC 



(III.67a) 
(III.67b) 
(III.67c) 



Here and in the following, we use certain combinations of both the magnetic 
three-vector potentials Aa{r) and of the electric potentials *"]4o(r): 



Ai 
An 
Am 



% + % 
% + % 
% + % . 



(IIL68a) 
(III.68b) 
(III.68c) 



A + As % 
A, + As 
A, + A, 

The reason for this is that these combinations transform inhomogeneously as 

Ai^ A[ = Ai + Vai (III.69a) 

(III.69b) 
(III.69c) 



All A[^ = All + Va2 
^iii All = An + Vas , 



provided the wave functions ipa{^) and the exchange potential B{f) un- 
dergo a homogeneous magnetic gauge transformation of the following kind 
(a = 1,2,3): 



B{f) =^ B\r) = e-*("2-"3)^(^) 
Bo{f) =^ B'^{r) = e-*("2-"=*)So(f) . 



(IIL70a) 
(111.70b) 
(111.70c) 



(For a general discussion of the RST gauge transformations see ref. [T2].) 
As a consequence of this gauge arrangement, the negative Pauli components 
^"'(/9_(r) (1111.67a|l - (1111.67c[) do inherit the same homogeneous transformation 
behaviour, i. e. 



(a), 



'V?-(rJ 



(a), 



'v?_'(f) 



-iaa . (a) 



'(/?_(f) , 



(111.71) 



so that one can expect that the emerging energy eigenvalue equations for the 
positive Pauli spinors (r) will also be found to be gauge covariant. 
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And indeed, substituting those negative Pauli spinors (r) (lIII.67a| - 
(lIII.67c| into the relativistic mass eigenvalue equations (lIII.40a| - (lIII.40c|) 
yields the corresponding non-relativistic linearized energy eigenvalue equa- 
tions for the positive Pauli spinors (r) in the following form 



2M, 

re 



- (V + iA,Y + he % + ^ (^i • ^) } = i^p(i) ■ 

(III.72a) 



2^ (V + ^in) -nc^%^^[H.rB) 



> 2M ^ ' ° 



2M, L 



(V-S) +2S- V 



(V = Ep(2) ■ (V 

(III.72b) 



2M, 



(V + ^AnY - he + ^ (#in • a) } ®v._ 

(v-i) +2i-v 



h"^ i * h^ 



2M, 



Ep(S) ■ • 
(III.72c) 



Here the (non-relativistic) Pauli energy eigenvalues -E'p(a) are defined in a 
nearby way in terms of the particle rest masses Mp, Mg and the mass eigen- 
values Ma through 



E, 



p(i) 

'P(2) 
'P(3) 



-(Mp + Mi)c2 
-(Me - M2)e^ 
-{M,-M;)e^ 



(III. 73a) 
(III.73b) 
(III.73c) 



Observe again the specific way in which the Pauli eigenvalue system 
flIII.72a[) - flIII.72c| l reflects the interactive structure of the three-particle sys- 
tem: Equation (IIII.72a)) for the first (different) particle says that this particle 
is subjected exclusively to the electromagnetic interactions with the other 
two (identical) particles, namely via the electric and magnetic potentials 
(!III-68a| . where the latter potential Ai refers to the interaction due 
to the translational (point particle) degrees of freedom. On the other hand, 
the presence of the magnetic field f^i (= V x Ai) in equation flIII.72a|l is due 
to the magnetic dipole character of the first particle. As a consequence of 
this ordinary interaction via the electromagnetic potentials '^o and Aj, the 



42 



first eigenvalue equation (lIII.72a[) is of tlie usual Hamiltonian form 



m 



Epil) ■ %+ 



where the Hamiltonian Hi is evidently given by 



Hi 



2M 



{y + iAiY + hc^A, + ^{Hi'a) 



(III.74) 



(III.75) 



Clearly, such a purely electromagnetic type of interaction is active also 
for the second and third particle, cf. equations (1III.72b[) - (]III.72c|) : but ad- 
ditionally these two {identical) particles do feel the exchange forces, where 
the latter subdivide again into those of the electric type (~ Bq) and those of 
the magnetic type (~ B). Consequently, the eigenvalue equations flIII.72b l- 
(IIII.72c|l for the two identical particles are not of that simple form as for the 
first particle flIII.74[) but rather do appear in the modified Hamiltonian form 



Hii%+ + hii%+ = E, 
Hiii%+ + hiu^W = 



P(2) ■ 

P(3) ■ %^ 



(HI.76a) 
(ni.76b) 



Here the ordinary Hamiltonians Hu^m are of the usual electromagnetic type 
fHirfSjl . i.e. 



ILIII 



■^(V + zi'„,„:)^-;icWiK + ^(#n,in-a) , (III.77) 



and furthermore the "exchange" Hamiltonians hu m can be read off immedi- 
ately from flIII.72b)) - (llll.72cp as 



hii 



h 



III 



2Me 
2M, 



(Y-a) -hcBo 
{Y-a) -hcBo 



2M, 



2M, 



[V'B) 

(v-i) 



2B-V 



2B'V 



(III.78a) 
(ni.78b) 



It is true, the Pauh eigenvalue system (]III.72a[) - (]III.72c|) is not strictly 
linear because the terms quadratic in the vector potentials ^i,ii,iii have been 
retained; but this is only in order to formally ensure the gauge covariance of 
that Pauli system where the exchange Hamiltonians flIII.78a|l - (IIII.78b) must 
be presumed to transform homogeneously, as is the case with the exchange 
vector potential B ( 1III.70b[) . It should be a matter of course that the strict 
gauge covariance can exist only in the original non-linear theory! 
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5. Pauli Energy Functional 



Once the non-relativistic approximation of the original relativistic mass 
eigenvalue equations flIII.40a[) - flIII.40c|l is now firmly established, one would 
like to go also the other way round by recovering just these non-relativistic 
equations as the variational equations due to the non-relativistic limit of the 
original mass functional M^tjc^ flIII.50[) . However, this non-relativistic Pauli 
energy functional (-Epj, say) can easily be obtained by simply looking for 
the non-relativistic limit forms of each of the constituents of M[t]C^. This 
procedure will also further elucidate the mathematical consistency of the RST 
formalism. From the physical point of view, it should be self-evident that the 
practical usefulness of such a non-relativistic functional -Epj will consist in the 
development of variational techniques for obtaining approximative solutions 
to the non-relativistic eigenvalue system (lIII.72a)l - (lIII.72c|) . see below. 



N on- Identical Particle (a = 1) 

As the first one of these limit forms consider the relativistic kinetic en- 
ergy Ikin, which is the sum of the individual single-particle contributions 
Tkin(a) (1111.461) - flIII.47[) . Eliminating here the negative Pauli spinor ^Ap_{f) 
from the kinetic energy 7kin{i) of the first particle (IIII.46| by means of the 
approximation (lIII.67a[) yields 



kin(l) 



in{l) 



2M„ 



(1) 



+ i 



^ /'dV(V{'?-(V + .Ai),a-Ai}(V. (111.79) 

P J 



But here one can resort to the operator identities 



a'{V + iAi) ={V + iAif -B'Hi (111.80a) 
|a- {V + iAi),a' Ai^ = V • li + 2v4i • {V + iAi) +ia'Hi (111.80b) 
(with V • i*! = 0) , 
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and this recasts the non-relativistic hmit i^ji^^ of Ikin(i) (1III.79| into its final 
form 

+ 4|-/dV(V('?-^i)^V-y /dVlrq,. (III.81) 

The non-relativistic approximation q^^ of the convection four-current gi^ fllll.54[) 
is easily found here as 

cii = -^{(V(V + ^iI)(V-[(V + ^lI)(V]'^V}, (111.82) 

namely by simply omitting the negative Pauli spinors ^^\p- for its spatial part 

gi- 
lt is interesting to remark here that the non- relativistic limit E^^^^^^ f llll.Sip 
of the kinetic energy 7kin{i) does not agree with the expected form i4;in(i) from 
conventional quantum mechanics, i. e. 

^in(i) = -^y"dV-'(V(V + ^A)^(V, (111.83) 

rather this does appear only as the first part of i^ijj(i) f!II1.81| ! Obviously, 
the relativistic kinetic energy 7kin(i) (1III.46|I does contain also the magnetic 
interaction energies of the polarization and convection type (-w second line 
of equation (1III.81[) ) . This is of relevance for the non-relativistic limit of the 
first particle's contribution M[i]C^ to the total mass functional M[t]c2 f llll.43p 

-M[i]c2 = Mpc^ ■ Z^i) + 2rkin(i) + 2Mf5"')c2 (III.84) 

because the magnetic interaction energies (on the right-hand side) must com- 
bine with the kinetic energies into a consistent non-relativistic approximation 
E[p] of the mass functional. For instance, consider the mass renormalization 
factor Z"^-^^ of the first particle which is given by equation flIII.44)l for a = 1. 

Eliminating from that equation the negative Pauli spinor ^\)_{f) by means 
of flIII.67a| with regard of the relativistic normalization condition (IIII.42)) 
yields 



Zfi) = 1-2 jd?r^V_%_ 
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and thus one finds by use of the operator identity (lIII.80a[) and the non- 
relativistic kinetic energy -£^in(i) (1III.83[) 

Therefore the first two contributions to the first mass functional M[i]C^ (1III.84[) 
combine to the following non-relativistic form 

Mpc^ ■ + 2rki„(i) Mpc^ + - he jd^fA • qi , (111.87) 

which evidently contains still the first particle's magnetic interaction energy, 
but now exclusively of the convection type (~ Ai) and no longer of the 
polarization type (~ Hi), as was originally the case with f lIII.81[) . 

Finally, for the non-relativistic limit of the mass functional M[i]C^ (IIII.84|) 
it becomes necessary to consider also the non-relativistic limit of the first 
particle's mass equivalent M^^^^c"^ flIII.52)) . whose splitting into a convection 
and a polarization part has already been specified in relativistic form through 
equation ( 1III.57|) . Turning here first to the convection part M^f^c^ 

Mgc^ = Jd'fqi.iA'^ + A'^) , (III.88) 

it is immediately obvious that this object splits up in a natural way into 
a sum of an electric time) part M^^^^c^ and a magnetic (— space) part 

mI^^'^c^ i.e. 

Mgc^ = V + Mf^'^^c^ , (III.89) 

with the electric part being given by the time components of the four-vectors 

as 

MgV = -^jd'r %{(% + (%) (111.90) 

and analogously for the magnetic part in terms of the corresponding space 
components 

Mfj-^c^ = ^ jd^Tq,'{A, + A,) . (III.91) 
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Clearly the non-relativistic limit of this magnetic convection part is easily 
written down by means of the non-relativistic convection current (1111. 82|) 
as 



he 



y"d¥qi-(A2 + i3) (111.92) 



and thus compensates for the magnetic convection energy in the non-relativistic 
limit ( 1111-81|) of the relativistic Tkin{i)- But also the electric convection 

energy Mfjj'^c^ (1X11901) receives a very simple approximation, namely by ob- 
serving that the time component of the convection current qi^ ( 1111. 54|) 



contains the time derivative of the first Dirac spinor V'i(r) (!ITL9) such that 



%{f) = 2 ( ^ - % ) MrlMr) -VArY'Wir) (111.93) 



2Mp 2Mp 



where for the non-relativistic limit one may put here Mi =^ —Mp. Thus the 
electric convection energy (llll.90|) becomes in the non-relativistic limit 



MfjV ^ m[^]V = ^ y"dV%(f)(V(O^V(0 , (ni.94) 

which together with its magnetic counterpart M|*j^j™'*c^ fixes then the non- 
relativistic limit of the total convection energy M^^fc^ (UTTsg) . 



Therefore, what finally remains to be determined (for the electromag- 



netic mass equivalent mIj^^c"^ (1111. 52|) ) is its polarization part mB^c^ (1111.58 



But this is a simple matter when one observes that the electric dipole den- 
sity ^"^Si (1111.59a|l of the first Dirac field is built up by the product of posi- 
tive and negative Pauli spinors ^\>±{r), whereas the magnetic dipole density 
'"'iSi (llll.59b|l contains also the product of the positive Pauli spinors ^^^ip-+{f) 



alone. If now all the terms with negative Pauli spinors are dropped for the 
non-relativistic limit, one arrives at the following limit forms of the Dirac 
dipole densities ( llll.59a)l - (1111.59b)l : 



^"^Si (III.95a) 

'-'Si ^ (^^^V , (ni.95b) 
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and this then fixes the non-relativistic hmit of the spin polarization energy 

Mf^jc^ furmii as 

M^/ Mgc^ = ^|dV(i?, + ^3)-(WV). (111.96) 

Summarizing, the electromagnetic mass equivalent of the first particle 
pil.52| becomes in the non-relativistic limit with regard of the present results 
mim^ . fUTLM)) and fUTLOG]) : 

+ ^|dVi?r(WV). (111.97) 

Of course, this result can be obtained also via the direct splitting of Mjy™-*c^ 

into a non-relativistic electric part m|^^jC^ and a non-relativistic magnetic part 

M[y^c2, cf. dnLH. But if this is now combined with the non-relativistic limit 
of the rest mass plus kinetic energy ( 1111.87)1 into the first particle's total mass 



functional (1111.84)1 . one finally arrives at the desired Pauli energy functional 
Epiiy. 

+ hc |dV%(V(V + ^ ld'fH,'{%la%^) . 

(III.98) 

And indeed, if this functional is further modified according to the method of 
Lagrangean multipliers as usual to £'p[i] 

-E^p[i] =i= -£^p[i] + Ap(i) Nd(i) , (111.99) 

with the non-relativistic limit Nd(i) of the relativistic normalization con- 
straint (IIII.42)) being easily determined as 



Nd(i) = jd^f%l%+-l = 0, (III. 100) 
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then the variational equation 



6E, 



m 



(III.lOl) 



with respect to the first Pauli spinor ^\)^ is found to just agree with the first 
particle's non-relativistic equation ( III. 72c^ ! Of course the first Lagrange 
multipher Ap(i) must be identified here with the first Pauli energy eigenvalue 



Epn) flIII.73a|) . i.e. 



A 



p(i) 



p(i) 



(III.102) 



It should also be a matter of course that the additivity of the relativistic 
mass functional M[t]C^ (IIII.41[) will transcribe to its non-relativistic limit, i. e. 



El 



[P] 



Ep[i\ + Ep\2\ + E] 



P[2] 



'P[3] ■ 



(III.103) 



But since the first Pauli functional £'p[i] (1III.99|) does depend exclusively 
upon the first Pauli spinor but not upon the other two spinors 
and ^\>+, the extremalization of i?[p] with respect to ^^\p+ is equivalent to the 
extremalization of merely the first Pauli functional Ep^i] (1III.99|I . Therefore it 
is not necessary for the variational deduction of the corresponding first non- 
relativistic equation flIII.72a| to determine the other two functionals E'ppj and 



61E1P] = SiEp[i] = 



(III.104) 



Identical Particles (a = 2, 3) 

However, such a variational independence cannot exist between the sec- 
ond and third Pauli functionals Ep[2] and E'pisj because the corresponding 
eigenvalue equations flIII.72b|) and flIII.72c)l display a direct exchange cou- 



pling between both Pauli spinors ^\)+ and ^%:>+. Therefore one expects that 
the desired Ep[2] will depend not only upon but also upon ^%:^+, and vice 
versa for -Eppj. But on the other hand, the general structure of both the 
original relativistic mass functionals (a = 2,3) 

M[a]c' = ZlyM,c' + 2T^,^(^,^ + 2^41X^0' -2Ml^fc' (III.105) 
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demonstrates that they formally differ from the first functional — M[i]C^ 
(1111.84)1 merely by the presence of their exchange parts M^^^^c^, i.e. 

(III. 106a) 

M^f^ = ^ Jd'rBo{f)ho{f)-^ jd?rB{f)'h{r) = M^/ - M^f\ 

(111.106b) 

while the first three contributions are quite analogous to the case of the first 
particle (a = l),cf. pil.84| . Therefore it is merely necessary to determine 
here explicitly the exchange contribution m'^^c^ for a = 2, 3. 

The "electric" part M^^c^ hereof describes the exchange interaction of the 
electric type and is given by 

Mg^c' = y Jd^f Bo{r)ho{r) (III. 107a) 

K]""' = Y /d'^"*^o(r)4(f) , (III.107b) 

where both contributions (a = 2, 3) have already been proven to be identical, 
cf. (1III.27[) . But here it is a very simple thing to determine the correspond- 
ing non-relativistic limit, namely by recalling the general definition of the 
exchange current /i^ ( III.24[) whose time and space components {ho; —h} read 
in terms of the Pauli spinors ^"'^ip±{f) 

h, = (3)^^ + (Ill.lOSa) 

h = (%la%+ + ^^^ifla^. . (Ill.lOSb) 

Thus the non-relativistic approximation Iiq of the time component ( ]III.108a| 
is immediately evident: 

ho ^ lio = , (III. 109) 

which then yields for the exchange mass equivalents of the electric type 
f llll. 107a)) - fllll. 107b)) 

M^lc^^ul^f^ = y jd^r Bo{f)%l{ffW{r) (Ill.llOa) 
Mj^^c^ ^ Mj^jc^ = y jd'fBo{ff\pl{rfW{f). (Ill.llOb) 



50 



But somewhat more intricate is the determination of the non-relativistic 
hmit of the exchange three-current h{f) (]III.108b[) : Substituting here the 
negative Pauh spinors ''^'(/:'_(f) from the approximations flIII.67b|l - flIII.67c|) 
one is led to the following splitting of the exchange current h 

h ^£ = b + i, (III.IU) 
with the (complex-valued) convection part b being given by 

b = ^ { (V(V)'^V - ^V(V(V) } . (111.112) 

Clearly this is just the non-relativistic linearized form of the exchange con- 
vection current flIII.64a[) and thus is the exchange counterpart of the elec- 
tromagnetic convection currents q„ (see, e.g., equation (]III.82)) for a = 1 or 
(]III.118a|) - (lIII.118b[) below for a = 2,3). Quite analogously to the electro- 
magnetic polarization currents Sa (e.g. equation (]III.61[) for a = 1), there 
emerges here the exchange polarization current w, as the curl of the exchange 
polarization density ^™'Z 

i = Vx("^)Z (111.113a) 

'-)Z = ^ (2U<?(=^W , (ni.ii3b) 

which is again the linearized non-relativistic version of fllll.64b[) - pil.65|) . 

Now as a consequence of that splitting of the exchange current h (1III.111[) 
into a convection and polarization part, the exchange mass equivalents ^j^fjc^ 
flIII.106a)) - flIII.106b[) do naturally split up also into a convection and a polar- 
ization part: 



M^[c' = — j&'r B{f)'h{r) M|c' = — J d'f B{f) 'h{f) 

+ 4|-/dVy-(WV) 

(III. 114a) 

M®c2 = y jd^rB{f)'h{f) M®c2 = y jd?fB{r)'h{r) 

(III.114b) 
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But the point with these non-relativistic hmits of the magnetic exchange mass 
equivalents is now that they just generate the required exchange terms for the 
non-relativistic energy eigenvalue equations flIII.72b[) - (IIII.72c[) . or flIII.76a l- 
flIII.76b|) . resp., i.e. exphcitly 



2] 



(111.115a) 



(III.115b) 



Now that the exchange terms are clarified, we are left with the problem of 
writing down the kinetic and electromagnetic parts of both functionals M]^a]C^i 
a = 2, 3 (1III.105[) . This problem, however, can easily be managed by simply 
referring to the analogous case of the first particle flIII.84|l . Namely, despite 

the presence of an additional exchange term (-^ B, B) in the approximate 

expressions (]III.67b|l - (lIII.67c|) for the negative Pauli spinors ^\p^ and ^^V-) 
the rest mass and kinetic energy of both identical particles (a = 2, 3) do 
combine in the linear approximation to a result quite analogous to the first 
case f lIII.87|) . i.e. 



MeC^ ■ Zf2) + 27J;in(2) 



'kiii(2) 



he / d f All * q2 



M^c' + £iin(3) - he /dV Am • qg 



(III.116a) 
(111.116b) 



where the non-relativistic kinetic energies £kin{a) ^oi both identical parti- 
cles (a = 2,3) are of course defined quite analogously as for the first par- 
ticle fUTLSSjl : 

^kin(2) = [d'f^\l{V + tliiY'^'W (III.117a) 



El 



'kin(3) 



2Me 



jA'f%i{v + iAiiy%^ 

|dV(^U(V + zl„i)'(V, 



and similarly for the non-relativistic convection currents (a 



ih 



q2 ^ 



qs 



2MeC 



ih 



2M,c 



%i (V + lAiii) (V - [(V + lAiii) % 



(2), 



(3), 



(III.117b) 

2,3) 

(111.118a) 
(III. 118b) 
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cf. (1III.82| for the analogous case with the first particle. (However, observe 
here that the Gordon decomposition for the identical particles a = 2, 3 reads 

Kf, = Qai^ + Sai,, (III.119) 

in contrast to the case ( 1111.53)1 with the first, different particle a = 1. See 
also below for the three- vector version fllll.l23|) hereof.) 

But also the electromagnetic contributions M^^^^^c^ of both identical par- 
ticles are nothing else than copies of the first case Mf5°'^c2 (IITLQT)! . However, 
it is very instructive to deduce this also alternatively via the original de- 
composition fllll.52)) . which refers to an electric (Mr^^c^) and a magnetic part 



(M^l^c^), rather than via the other decomposition fllll.57| into a convection 
fj^j^c^) and a polarization part {M^^f-^c 



part (Mr^j^c^) and a polarization part (Mr^^jc^). Here the electric contributions 



Mj^^jC^ appear especially simple because one merely has to approximate the 
Dirac densities ^°'^ko{f) by the first term with the contraction of the positive 
Pauli spinors '"V+IO^ i-^- 

(% = ('^v+V + 'V'V = '"V+V, (111.120) 

and this immediately yields the non-relativistic limits of the electric mass 
equivalents (1111. 48a[) as 



(111.121a) 



~ jd?r^'%{r).%l{r)%^{r). 



(e) ^ 



(111.121b) 



However, the analogous treatment of the magnetic mass equivalents m'^^c^ 
fllll.48b[) is somewhat more complicated. First observe here that for their 



linearized form the negative Pauli spinors '"'99_(r) (1111.67b|l - (llll.67c|) may be 
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further approximated as follows 



-^J' (V + zAiO^V (III.122a) 



(3), 



ih 



2M,c 



a ' {V + iAui)%+ . (III. 122b) 



Next by use of these approximations, the Dirac currents ka{r) are split up 
into the sum of a convection part and a polarization part 

ka = V<?'V + *"V|'?'V ^ ia = ^a + Sa, (HI. 123) 

where the convection parts are given by flIII.118a[) - flIII.118b|l and similarly 
the polarization currents are found again to be of the analogous form as for 
the first particle ( ]III.61|l . i.e. 

s„ = V X ^'^^Sa (III. 124a) 

'"'Sa = ^^M^^^W (III. 124b) 

(a = 2,3) . 

Consequently, if this splitting of the Dirac currents ka{f) (1III.123|) is sub- 
stituted into the magnetic mass equivalents M^^^c^ ( ]III.48bll there occurs a 
corresponding decomposition into two parts: 



'[2]- - 2 



(III.125a) 



m(?c2 = /d¥4-An^Mf?c2 



[3] - - 2 J ' ' 

(III. 125b) 

Clearly it is not surprising that these magnetic mass equivalents of both 
identical particles appear very similar to the case of the first particle, cf. ( ]III.97)l . 
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which then apphes also to the combined electric and magnetic objects 



(III. 126a) 



2 



4M, 



dVi?ni-CWW)- 
(III.126b) 



This agreement of the non-relativistic limits of all three electromagnetic ob- 
jects m'^^c^ (a = 1,2, 3) validates the former assertion that this limit is 

independent of the special decomposition of m'^c^^ i.e. either the decom- 
position with respect to an electric and a magnetic part (as shown by equa- 
tion flIII.52| ) or with respect to a convection and a polarization part (as 
shown by equation (]III.57|) ). 

But the crucial point with these electromagnetic mass equivalents is now 
again that, through their combination with the kinetic and rest mass terms 
(IIII.116a)l -( ]III.117b|) . the convection parts (~ q^) do cancel so that the non- 
relativistic Pauli functionals of both identical particles (a = 2, 3) as the 
non-relativistic limits of the mass functionals Mya]C^ (1111.105)1 

M^a]C^-M,^ ^ Ep[a] (III. 127) 

(a = 2,3) 
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are emerging in a relatively simple shape, namely 

En2] = /dV(^U(V + - he jSf^X . 

+ ^ Jd'fHu' - 2M^f (III.128a) 

Evidently, these Pauli functionals for the identical particles are of the same 
form as that one for the different particle pil.98|) , apart from the additionally 
emerging exchange mass equivalents M^^^c^ ( ]IlI.106a|) - (1III.106b)l which can 
of course occur only for identical particles. But with the present Pauli energy 
functionals fllll.l28a|) - fllll.l28b[) our final goal has now actually been attained; 
namely, the variational equations 

5Ep[a] = (111.129a) 

Epla] = Ep[a] + Ap(,) ND(a) (111.129b) 

(a = 2,3) 



turn out to be identical to just those non-relativistic Pauli equations f lIII.72b l 



(lIII.72c| . Here the Lagrangean multipliers Ap(a) are to be identified again with 
the Pauh energy eigenvalues (]III.73b|l - (lIII.73c[) : 

Ap(a) = -^P(a) , (111.130) 

and the non-relativistic approximation to the normalization constraint (1III.42[) 
is of course given by 

Nd(.) = ydV'^Vi'V-l = 0. (111.131) 

Summarizing the results, the relativistic mass eigenvalue equations (lIII.40a[) 
(IIII.40c| actually do admit the interpretation of being just the variational 
equations due to the mass functional M[t]C^ fllll.50[) . which is the sum of the 
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three individual mass functionals M[a]C^ (a = 1, 2, 3). This sum is decoupled 
in the sense that the ath eigenvalue equation of the system (lIIL40a| - (lIII.40c|) 
is the variational equation (with respect to the ath Pauli spinor *"V±) of the 
ath mass functional M[a]C^. The non-relativistic limit of this mass eigen- 
value system is given by equations flIII.72a[) - flIII.72c|l . which can be obtained 
in two different ways: (i) either one eliminates the negative Pauli spinors 
''^'(yC„(r) directly from the relativistic eigenvalue equations by means of the 
approximations (]III.67a|) - (lIII.67c[) . or (ii) one first eliminates the negative 
Pauli spinors from the relativistic mass functional M[t]C^ in order to obtain 
its non-relativistic limit, i.e. the Pauli energy functional E[p] ( 1III.103[) . and 
afterwards one looks for the variational equations due to this non-relativistic 
functional E[p]- Thus the logical scheme sketched at the beginning of sub- 
section 137] of the present section has actually been validated. 

But though this coincidence of the results (obtained along different routes) 
is very satisfying, one nevertheless cannot be contented with this situation 
since it is obvious that the gauge field equations (e. g. the Poisson equations 
f lIlL3) - flITL4) l) cannot be deduced via a variational procedure from that rel- 
ativistic mass functional Af[T]C^ nor from its non-relativistic approximation 
E[-p]\ The reason is here that both the mass functional M[t]C^ and its non- 
relativistic approximation E'p] do not contain any derivative of the gauge 
fields! Therefore, if we wish to possess a (relativistic or non-relativistic) en- 
ergy functional, whose variational equations are required to embrace both the 
matter and gauge field equations, then we have to resort now to a new con- 
struction: the RST energy functional (i?[T])- However, the mass and Pauli 
functionals discussed so far will thereby not become obsolete, but on the 
contrary will be of valuable help for constructing the desired RST functional 



IV RST Energy Functional 

After it has become sufficiently clear now that the mass functional M[t]C^ 
cannot directly serve as the desired energy functional i?[T], one is motivated 
to turn back to the total RST energy (1II.50| . Indeed, this object contains 
both the matter energy E'd (III.52a| and the gauge field energy Eq (111.52b|l . 
which is a necessary condition for the variational deduction of both the mat- 
ter plus gauge field equations from the corresonding functional E[t]. Never- 
theless, this total energy functional i?[T] cannot be immediately adopted as 



57 



the desired goal. The reason is that the matter and gauge fields are coupled 
in such a way that the mass eigenvalue equations plus Poisson equations do 
not appear as the corresponding variational equations! On the other hand, 
the definition of the field energy Et flll.50|l as the integral over the energy 
density *'^'3oo(r) appears to be very natural and absolutely consistent with 
all that one knows about the logical structure of the successful field theories 
in theoretical physics. Consequently, when one wishes to modify the original 
RST functional E[t] in order to meet with the complete set of field equations, 
then one must try to admit only those modifications which do preserve the 
numerical value of the original E^. A special class of such modifications con- 
sists obviously in adding to the original functional i?[T] some constraint which 
is automatically obeyed by all solutions of the RST eigenvalue problem. In- 
deed, we will readily demonstrate that constraints of this kind are provided 
by just those Poisson identities, cf. equations fllll.35|) - fllll.36|l . but not by 
the exchange identities (llll.33|l - fllll.34|l . so that the inclusion of the exchange 
forces requires a slight modification of the original proposal -Et flll.50[) . 



1. Physical Meaning of the Poisson Identities 

The role played by the Poisson identities becomes elucidated by a closer 
inspection of the splitting of the total field energy E^t] (111.51| . Its matter part 
£^[D] (111.52a| has already been shown to be the sum of three single-particle 
contributions Ej)[a] (a = 1,2,3) 

= £'d[i] + -E'd[2] + E^s\ , (IV.l) 

which then adds up to the total mass functional M[t]C^ minus the mass 
equivalents of the electric and exchange type [H [H] 

E[j)] = M[t]c2-2(M(V-M^)c2) . (IV.2) 

Recall here that the total mass functional has already been specified by equa- 
tion flIII.43[) so that the matter energy flIV.2p appears now in the following 
(relativistic) form 

+ 2rkin - 2 (m("^)c2 - Mfe^c') . (IV.3) 

This is a very plausible result because it says that the energy of matter 
consists of the (renormalized) rest mass energies of the three particles (first 
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line on the right-hand side of (1IV.3[) ) plus the proper kinetic energy of motion 
27kin minus the magnetic mass equivalents of the electromagnetic (~ M^™)) 
and exchange (~ M^^) type. Observe here that the electric mass equivalents 
(~ M^^'^^) do not contribute to the matter energy Ed (IIV.3|I . quite in contrast 
to the mass equivalents of the magnetic type (~ M^™'^)) which represent 
here the field-theoretic analogue of the well-known "minimal substitution" 
(p — > p + hA) for the coupling of a spinless point particle of momentum p 
^V) to a three-vector potential A, see the kinetic energies -Ekin(a) (1III.83|) 
and kll.ll7a)l - flIII.117b)l . 

But concerning the Poisson identities, they are rather connected with the 
gauge field contribution Eq flll.52b|) of E^. This is better reahzed by splitting 
up Eg into its electromagnetic and exchange parts, i. e. 



Er 



En — Er 



(IV.4) 



where the energy content of the real modes (~ ^4"^) is built up again by 
the electric (e) and magnetic (m) contributions [8J 
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and similarly for the complex modes 

Ec 
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(IV.5a) 
(IV.5b) 

(IV.5c) 

(IV.6a) 
(IV.6b) 

(IV.6c) 



Observe here the important fact that the exchange energy Eq flIV.6a|) - flIV.6c[) 
is always positive [Eq > 0) but enters the gauge field energy Eq f lIV.4[) with 
a negative sign. This implies that it becomes energetically more favourable 
to undergo the exchange interactions which, on the other hand, is possible 
only for identical particles if the basic conservation laws are to be respected, 
see below! 
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The relationship between the Poisson identities and the various con- 
tributions (lIV.5a| - (lIV.6c| becomes now immediately evident through sim- 
ply recalling the definitions of the field strengths in terms of the poten- 
tials, cf. flllI.15a|) - flIII.15d|l . For instance, consider the electric Poisson iden- 



tity (II11.35|) and find that this identity expresses nothing else than just the 
identity of the electric gauge field energy E§ (IIV.5b)l and the electric mass 
equivalent M^^^c^ (llll.48a)l . i.e. 

iV|) = Ef-M^V = 0. (IV. 7) 

Or in other words, the electrostatic gauge field energy equals the elec- 
trostatic interaction energy M^^^c^ of the matter field modes with the gauge 
field modes. In a similar way, it is easily seen that the physical content 
of the magnetostatic Poisson identity (1III.36| consists in the identity of the 
magnetic gauge field energy E^^^ (dVlcJ and the magnetic interaction energy 
M^(? (lIII.48b)l of the matter and gauge field modes, i. e. 



AtH = _ M^^V = . (IV.8) 

However, the situation with the analogous exchange identities is some- 
what more intricate. Namely, the field equations flIII.10a|) - flIII.10b)) for the 



exchange potentials Bq, B are not of the simple Poisson form but rather 
contain a "mass term" (~ a^) and therefore do rather resemble the static 
form of the Klein-Gordon equations. When this fact is combined with the 
former integral relations (]III.31[) one arrives at the exchange counterparts of 
the electromagnetic identities (1IV.7[) - (1IV.8[) in the following form: 

ATf) = eI") + ■ \- [d^fB^^B^ - M^%^ = (IV.9a) 

iV® = + • ^ /"dV B^B^ - Mfe)c2 = . (IV.Qb) 

This result says that the energy content of the exchange field modes (of the 
electric (h) and magnetic (g) types) differs from the corresponding exchange 
mass equivalents M'^^c^ (1111.27)1 and M^^c^ (1111.28)1 by the spatial integral of 
the Lorentz invariant B^^B^. Thus it is only for states with infinite exchange 

length om flIII.7)l that the exchange energies e'^ and e'^ (lIVJb]|-(IIVJc]| can 
coincide with the corresponding mass equivalents M^^c^ and M^^c^. 
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This non-identity (for finite om) of the exchange energies and their mass 
equivalents will readily force us to modify the RST field energy (III.50| in 
order to obtain a physically reasonable energy functional i?[T]- 

2. Relativistic Variational Principle 

The construction of the desired relativistic energy functional £'[t] has 
already been described in great detail for a system of two different particles 
(i.e. positronium, ref. |9j). Therefore it may be sufficient here to briefiy 
comment on the analogous construction for the present three-particle system 
(of two identical and one different particle). Namely, the principal idea refers 
again to the modification of the original RST energy E^ (III.50[) by means 
of the aforementioned Poisson and exchange identities so that the numerical 
value of the total energy is preserved (as far as possible). 

In the first step, one puts together the total energy i?[T] (III.51|) - flII.52b[) 
from the matter E'p] flIV.2| and the gauge field energy E[q] flIV.4)) - flIV.6c| in 
order to find the first tentative proposal E^!^^ in the following form: 

= M[t]c2-2(m(V-M^)c2) +4 + 4°^)-4)-4. (IV. 10) 

Next, the mass equivalents of the electric type are eliminated from this first 
proposal just by means of both identities (1IV.7|) and (lIV.9a[) of the electric 
type, which then yields the second tentative form as 

= M[^]c'-E§ + E^^^+eS^-E^S^+^.^ [d'fB^B,. (IV.ll) 

From this intermediate result a third proposal is obtained now by substituting 
here the mass functional M[t]C^ from equation flIII.43|l . which brings the total 
RST energy functional -E[t] to its final form 

+ E§ - - 4) + 4 + ■ ^ [d^r-B^B, . (IV.12) 

Observe here that the numerical value of the corresponding functional 
upon the solutions of the RST eigenvalue problem will exactly repro- 
duce the RST energy i?T in its original form ( jll.SOjl ! The reason is that all 
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those manipulations leading us from that initial form (1II.50| to the present 
E'ppj (IIV.12| are exactly valid for the RST solutions. Furthermore, the latter 
result fllV.12)) appears to be very plausible because the first line refers to the 
matter subsystem and thus specifies its rest mass and kinetic energy while 
the second line obviously consists in the energy of the gauge field subsystem. 

But now comes a subtle point which just refers to this numerical equiv- 
alence of the modified functional E^!^-^ fllV.12| and its original predecessor 
E'Y flll.50|l . Namely, a very natural expectation is that the gauge field sub- 



system should contribute its part to the total energy E't via the field strengths 
(i.e. E,H;X,Y), but not via the potentials (i.e. '"Mq, A^; Sq, -B). But ob- 
viously, this expectation is violated by the last term on the right-hand side 
of (]IV.12[) . which is built up exclusively by the exchange potentials B^. On 
the other hand, the presence of that questionable term is indispensable from 
the mathematical viewpoint because it is a consequence of the rigorous use of 
the Poisson and exchange identities, which ensure the numerical equivalence 
of Et (III.50|) and E^^ (IIV.12|I . Thus the mathematical rigour runs into con- 



flict with the physical intuition. In order to escape from this dilemma, one 
needs some convincing argument which settles the conflict either in favour 
of physics or in favour of mathematics. Here it seems wise to take the prag- 
matic viewpoint; and this means that one lets settle the question through the 
circumstance wether the stationary gauge field equations (]III.3|l - (]III.4)l and 



(Illl.lOap - dlll.lOb)) can be deduced from the functional flIV.12| . either 



with the presence or the omission of the questionable potential term. An- 
ticipating the result, physics gets the better of both and the potential term 
must be omitted, i. e. we adopt the desired energy functional {E^,^j, say) to 
be of the following form 

Ef^\ = ) ■ MpC^ + Zfs) ■ MeC' + 2^3) ■ MeC^ + 2Tkin + 4'^ - 4"^ - 4^ + 4^ ' 

(1V.13) 

otherwise the exchange equations (lIII.10a|) - (lIII.10b[) could not be identified 
as the variational equations due to the corresponding energy functional. 

Nevertheless, the present result ( ]IV.13|) . if interpreted as a functional 
i^p^-* acting over the configuration space of matter and gauge fields {^'^V±(0) 
A"-^{f), B^{f)}, will not yet (but almost) reproduce the RST eigenvalue 
system. Recall here that the latter system consists of the Poisson equa- 
tions fllll.3| - fllll.4[) for the real gauge field modes A"-^, the exchange equa- 
tions flIII.10a[) - flIII.10b|) for the complex modes B^, and also of the mass 
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eigenvalue equations (lIII.40a| - (lIII.40c[) for the Pauli spinors ^"'^ip±{f). How- 
ever, in order that the members of this coupled system of matter and gauge 
fields do actually emerge as the variational equations of an RST energy func- 



tional {Em , say) one merely has to complement the present result i?™ flIV.13[) 



by certain constraints via the method of Lagrangean multipliers. Clearly, the 
first one of these constraints must refer to the relativistic normalization con- 
dition flIII.42|) . But also the exchange identities flIII.33[) - flIII.34[) and Poisson 
identities f lIII.35[) - (]III.36[) must be adopted as constraints for the variational 



procedure. Thus one ultimately arrives at the desired energy functional i?[T] 
in the following form: 



a=l 



(IV. 14) 

And indeed, it is now a very instructive and enlightening exercise to 
identify the whole set of RST field equations as just the variational equations 
due to that functional E^^y. 

1. the relativistic mass eigenvalue equations (lIII.40a|) - (1111.40c)l appear as 
the variational equations with respect to the Pauli spinors *"'(/?-j-(r) (or 
their Hermitean conjugates *"V±(0) resp.) 



dEi 



[T] 



. (IV.15) 



Here the Lagrangean multipliers AD{a) are due to the normalization 
constraints and are to be identified again with the mass eigenvalues 
Ma fllll.9| , quite analogously to the case of the mass functional 0III.51a| - 
f llll.51c|l : 

Ad(i) = Mic^ (IV. 16a) 

Ad(2) = -Msc' (lV.16b) 
Ad(3) = -Msc' . (IV. 16c) 



2. the electromagnetic Poisson equations (1III.3| - (1III.4| turn out as the 
variational equations with respect to the gauge potentials '"Mq, Aa- For 
the link between the Maxwellian charge and current densities ja 
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and the Dirac densities ^'^^kQ, ka see equations (lIL40a[) -( lII.40e[) . The 
Lagrangean multiphers due to the Poisson constraints must adopt the 
following values: 

= -2 (IV. 17a) 

Ag^^ = 2. (IV. 17b) 



3. the exchange equation ( ]III.10a|) of the electric type is found to play the 



role of the variational equation of i?[T] with respect to the variation of 
the exchange potential Bq. Here the Lagrangean multipliers Xq and 
Aq due to the exchange constraints (1III.33[) - (1III.34[) must be chosen in 
the following way: 

Ag^=2 (IV. 18a) 

A|^=-2, (IV. 18b) 

which completes the set of Lagrangean multipliers. 



4. the exchange equation (IIILlObj) of the magnetic type is found to appear 



as the variational equation of £'[t] with respect to the variation of the 
magnetic exchange field B. Both exchange equations (]III.10a[) - (]III.10b|) 
represent the crucial point for the preceding discussion of the potential 
term (~ B^B^) in the rejected energy functional E^^^ (H\n2|l . Indeed, 
it is easy to see that both the electromagnetic Poisson equations (1III.3 I 



(1111.4)1 and the mass eigenvalue equations (lIII.40a| - (lIII.40c[) are insen- 
sitive to the presence or omission of that questionable potential term. 
But surely it is rather satisfying that one has to omit this term on two 
mutually supporting grounds: namely, (i) physical plausibility and (ii) 
variational deducibility of the exchange equations flIII.10a)) - flIII.10b|l . 
However unfortunately, the original definition flll.50[) of the RST en- 
ergy i^T must be abandoned to a certain extent under this compulsion 
of making a decision in favour of physical plausibility. 

Observe here that the energy functional i?[T] (IIV.14|I consists of two rather 
different contributions, namely of the physical contribution E^^^^ f lIV.13| . 
which is the sum of all physical energies (i. e. rest mass, kinetic, electro- 
magnetic, exchange energy) and of the constraints (i. e. wave function nor- 
malization, Poisson and exchange identities). For the practical applications. 
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where one may resort to appropriate trial functions for the RST fields (wave 
functions and gauge potentials), one will therefore select the trial functions 
in such a way that all the constraints are obeyed automatically so that it 
is sufficient to extremalize the physical energy Efl^l f HVl3) rather than the 

constrained energy E[t] flIV.14|) . 



3. Non-Relativistic Principle 

For the concrete applications, e. g. in the field of atomic and molecular 
physics, one can frequently be satisfied with the non-relativistic approxima- 
tion of the desired predictions. Thus if the RST principle of minimal energy 
is selected as the preferred variational method, one will be interested also in 
its non-relativistic limit. Clearly, one expects here that the non-relativistic 
form of the variational technique can be handled more easily than its rel- 
ativistic version. In this sense one wishes to see now the non-relativistic 
approximation of the present energy functional E'fT] flIV.14[) . But naturally 
for this purpose it is merely necessary to look for the non-relativistic forms 
of the individual contributions of E[t] ■ 

Turning first to the physical contribution iJ^fj^ (1IV.13|) of £^[t] (lIV-141) . one 
recalls that the particle rest masses {MpC^ and MgC^) are to be neglected 
for the non-relativistic approach and thus the matter contribution to E^^^ 
(first line on the right-hand side of (1IV.12|) ) is simply replaced by the non- 
relativistic kinetic energy ii^in as the sum of the three single-particle energies 

-Ekin{a) 

-Ekin =!= £kin{l) + ^in{2) + -Ekin(3) (IV. 19) 



as defined by equations flIII.83[) and flIII.117a| - flIII.117b|l . Next, consider 



the energy content of the gauge field subsystem, which is given by the sec- 
ond line on the right-hand side of equation (1IV.12|I . Since the gauge field 



energy Eq ( ]IV.4[) is defined exclusively in terms of the field strengths, see 



equations (lIV.5a| - (lIV.6c| . this energy contribution remains formally invari- 
ant when passing over to the non-relativistic limit since this limit concerns 
essentially the motion of matter, but not directly the mechanism of genera- 
tion of the gauge fields (Maxwell equations). Therefore the non-relativistic 
approximations do refer preferably to the matter subsystem, which (apart 
from the kinetic energies) enters the energy functional i?[T] flIV.14| via the 
constraints of wave function normalization ND(^a) flIII.42| . Poisson identities 
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N^^ and N^^ flIII.35)l - flIII.36)l . and exchange identities N^^ and A^® f lIlL33 l- 
(IIII.34|l . Here the non-relativistic form ND(a) of the original wave function 
normahzations N^^a) has aheady been specified through equations flIII.131[) 
so that we are left now with the task to present also the non-relativistic shape 
of the Poisson and exchange identities. 

Since, however, these identities are built up by the gauge field energies 
{E^^, E^^\ and their mass equivalents {M^V, MHc2, M'^^h'^, M^c^} 

see equations (lIV.7[) - (lIV.9b|l . it becomes merely necessary to look for the non- 
relativistic forms of these mass equivalents: 



(e) 



G 
(m) 



(e) 



'G 
(m) 
'G 



E, 



R 







G 



he 



G 



dV 



(v4)-(v5o) - — 4^0 



M^'>e 



(g) 



(g) 



^ 47raQ 



2 /d¥5'^5^ - M®c2 = 



(IV.20a) 
(IV.20b) 

■2 = 
(IV.20c) 
(IV.20d) 



Here the desired non-relativistic expressions have already been specified by 
equations (IIII.121ajl -( !riI.121b)l and (IIII.125ajl -( !riI.125b)l for the electromag- 
netic mass equivalents Mj^^jC^ and Mj|^''c^; and similarly the exchange mass 

equivalents m||^]C^ and M®c^ are given by equations (IIIL110a|l - flIII.110b[) 
and flHI.114a| - flIII.114bL All these results can now be put together for 
the non-relativistic version (E[t], say) of the relativistic energy functional 
i?[T] (lIV.14p . which then ultimately appears in the following form: 



E 



[T] 



+ he I d f 

3 



+ Ap(,) ■ Now + Ag- n|)+ Ag). N^) + A® ■ n| 



(IV.21) 



o=l 



Surely this is a plausible result from the physical point of view because 
it says that the non-relativistic energy Et of any RST field configuration is 
composed of essentially three constituents: (i) the first line represents the 
sum of the kinetic energy £kin of matter plus all four kinds of the gauge 
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field energy, (ii) the second line is the interaction energy of the individual 
magnetic dipole densities '"iSa and the magnetic fields Hj u^m of the other 
particles no magnetic self-interactions); and finally, (iii) the third line 
is the collection of constraints. Observe here that the magnetic constraint 
n|"^ (]IV.20bjl does not explicitly emerge in the third (i. e. constraint) line 



because it has been eliminated in favour of the dipole interaction term (second 
line). As a consequence, the magnetic field energy appears now with a 
positive sign in the non-relativistic functional E[t] fllV.21)) . in contrast to the 
relativistic functional flIV.13| . 

Now it is again a nice consistency check to deduce the non-relativistic 
RST equations from that functional E[t] ( ]IV.21)I : 



the non-relativistic energy eigenvalue equations flIII.72a|l - flIII.72c[) are 



obtained through the extremalization of E[t] with respect to the posi- 
tive Pauli spinors (^^ = 1; 2, 3); 

the electric Poisson equations (]III.3|) are obtained through extremal- 
ization of E[T] with respect to the electric gauge potentials '"Mq, where 
merely the relativistic charge densities '"^o(0 (or '"'^o? resp.) are to be 
replaced by their non-relativistic forms ^°'\o{r), i.e. concretely 

(% = ("V+V + *V*V ^ = >|'V ; (IV.22) 



the magnetic Poisson equations (1111.4)1 emerge through extremalization 



of E[T] with respect to the magnetic gauge potentials Aa, where the 
relativistic current densities jaif) (or kai resp.) have merely to be 
replaced by their non-relativistic forms ka, see equation flIII.123[) for 
a = 2, 3 ( and similarly for a = 1); 



the electric exchange equations ( IIII.10a[) arise through the extremal- 



ization of E[T] with respect to the electric exchange potential Bq (or 
i?o, resp.), where the exchange density /lo(^) is to be replaced by its 
non-relativistic version ]ho(r), see equation fllll.l09[) : 



5. the magnetic exchange equation (IIII.10b|) is obtained through extremal- 
ization of E[T] with respect to the magnetic exchange potential B (or 

B, resp.). Clearly, the exchange current h is to be replaced here again 
by its non-relativistic form Ih 
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Summarizing, one actually has achieved the original goal, namely to con- 
struct a relativistic energy functional E[t] (together with its non-relativistic 
approximation E[t]) from the intrinsic logic of RST, such that the solutions 
of the coupled set of (non-)relativistic RST field equations do just represent 
the stationary points of those functionals! 



V Fermionic and Bosonic States 

Naturally, the establishment of a new variational principle must provoke 
intense endeavours in order to either falsificate the new theoretical structure 
or to validate it by testing it for various applications and thus elaborating 
its logical consistency and practical usefulness. Especially one first wishes 
to become convinced that the predictions of the new formalism do coincide 
with those of the conventional theory for some of the standard situations 
where the results are already well known and in perfect agreement with the 
observations. In this sense, we select the hydrogen state 2p3 for a comparison 
of the conventional and RST treatment. This one-particle example can be 
solved exactly in both the conventional theory [18] and in RST; and it should 
not come as a surprise that the results of both approaches do exactly coincide. 
For instance, the energy (£tonv) of the quantum state is identified with the 
mass eigenvalue (M2C^), see equation flV.35)) below. On the other hand, 
the proposed RST energy i^[T] flIV.14[) is composed of various contributions 



due to the matter and gauge fields to be complemented by the constraints. 
But despite this different physical status of both competing quantities .Econv 
and EiT], their values upon the selected state 2p3 do exactly coincide, see 

equation flV.53[) below. Moreover, the RST functional ^'[t] meets also with 
the plausible expectation that for the non-relativistic limit there should occur 
some kind of spin degeneracy. The reason for this is that the spin-orbit 
coupling must disappear in the lowest-order approximation where the spin 
effect upon the energy levels is neglected. Subsequently it is demonstrated 
that the RST formalism does correctly account for this non-relativistic spin 
degeneracy. 

But clearly, such a coincidence of the conventional and RST predictions 
for the one-particle systems can not yet ensure the general confidence in the 
claimed RST principle of minimal energy; but rather one would like to see 
this principle working well also in the field of the few-particle systems (or even 
many-particle systems). Therefore, in order to present a less trivial situa- 
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tion, we consider positronium as a typical but sufficiently simple few-particle 
problem. More concretely, we demonstrate the working of the principle of 
minimal energy as an approximative technique for calculating the energy of 
the bound few-particle systems (albeit only in the non-relativistic regime). 
The treatment of such a two-body problem is surely a more rigorous test 
of the principle of minimal energy since in RST the two positronium con- 
stituents (i.e. electron and positron) must occupy bosonic states [TT]; and 
if any physical relevance is claimed for that principle, it must be shown to 
produce meaningful results also for such exotic quantum states with non- 
zero boson number. However, in order to not plunge here too deeply in the 
tedious calculations, we are satisfied with a treatment of the state 2^Pi in 
the non-relativistic limit to be combined with the spherically symmetric ap- 
proximation. It turns out that already the choice of a simple trial function 
with only two variational parameters (3 and u (see equation flV.125|l below) is 
sufficient in order to predict the (non-relativistic) energy of the state up 
to an accuracy of 10% (see the figure on p. 11071) . This is the same magnitude 
of deviation as for the analogous groundstate predictions in ref. [H] and thus 
hints on the hypothesis that the non-relativistic RST predictions may come 
very close to their conventional counterparts (or even coincide with them) if 
more subtle approximation techniques were applied. 



1. Selection of Spinor Basis 



The Pauli spinors ^°'^f± (r) flIII.39|l are elements of a two-dimensional com- 



plex linear vector space which may be conceived as a representation space of 
the orthogonal group 5*0(3) (or its covering group SU{2), resp.). Accord- 
ing to the specific transformation law under this group, those elements are 
called "spinors". Any such spinor can be described in terms of its compo- 
nents relative to some spinor basis whose selection is arbitrary on principle, 
but an intelligent choice mostly facilitates the managing of the concrete prob- 
lems. Since the Dirac wave functions ipa (1III.39|) appear as the direct sum 
of two Pauli spinors '"V±5 one can select a separate 2-basis for any Pauli 

2-space. Especially for the description of the present groundstate situation 

i i i _i 

one ffist starts with the C basis {Co''^,Co' ^} for the "positive" Pauh com- 

11 1 _i 

ponent and its concomitant {Ci '^iCi ' ^} for the "negative" component, see 
ref. [11]. Each of these two basis systems consists of eigenspinors for total 
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angular momentum J {= L + S) such that 

j^cr = mhcr (v.ib) 

PQ'"^ = hH{l + l)Q''^ (V.lc) 
|2 = ^(s + 1) C/-'" . (V.ld) 

For the positronium groundstate to be treated here one adopts the quan- 
tum numbers j = |,/ = 0, m = ±| for the positive Pauh component and 
j = |, / = 1, m = ±1 for the negative Pauh component (s = | in any case). 
However, we do not use directly this basis system ( lV.la|) -( IV.ld)l but we equip 
it with a further phase degree of freedom by using the "uj basis" which is 
related to the above mentioned % basis" through [11] 



(+) 



1 1 



Cr (V.2a) 



^{-) _ g*60 . ^5- 2 (Y 2b) 



1 1 



,{+) 



Here 5 is a real constant to be fixed later on by means of additional require- 
ments. This constant b may (or may not) be different for different particles 
(i. e. b ^ bp for the positively charged particle (a = 1) and b ^ be for the 
electron (a = 2)). The new basis elements are still eigenspinors of the z 
component (J^) of angular momentum, i.e. 

J,4+) = -(S-^);i-4+) (V.3a) 
J,4-) = (V.3b) 
j^4+) = _(5_l)/^.^(+) (v.3c) 

lJ-^ = {b-^)h-u[-\ (V.3d) 



1 1 



■ cr (v.2c) 



1 



and thus coincide with the old basis system ( ]V.la|) -( lV.ld[) for b = 0. But 



for arbitrary b one has a new situation with a non-unique basis system over 
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three-space, i. e. 

4J(0 + 27r) = e^^2-«4J(0), (V.4) 

where {r, 'd, (p} are the spherical polar coordinates as usual. Especially for 
§ = I one has vanishing eigenvalues of angular momentum 

i4? = (V.5) 

(« = 5). 

This basis system (§ = |) is double- valued over the two-sphere S"^ and has 
been used for the treatment of the positronium groundstate in the preceding 
papers [9l HOl [11]; but for the sake of generality we leave the constant b 
unspecified for the time being. 

A first hint on how to fix the value of 5 arises now from the specific shape 
of the Dirac three-currents ka flIII.123| which decompose in spherical polar 
coordinates as 

K = <%-e. + <%-e^ + %-e^. (V.6) 

Decomposing here also the (relativistic) Pauli spinors '"'(/2j_(r) 0III.39[) with 
respect to the selected u basis flV.2a| - flV.2d|) as (a = 1, 2) 

(V(r) = <'^'7^+ ■ 4+^ + ■ (V.7a) 
(V (r ) = -2 1 <'^'7^_ • uj[^^ + ^"^S^ ■ uj[-^ } (V.7b) 

lets appear the Dirac (relativistic) densities ''^'A;o(0 (1III.120|) in the following 
form 

<%(f) = . (V.8) 

Evidently, these densities do not yet provide an immediate handle for fixing 
the parameter 5. 

This situation changes now when one considers also the Dirac currents 
ka ( IV. 6)1 ■ which by their very definitions fill. 43a|l - f ill. 43c[) are always real- 
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valued objects: 

Wfc^ = — (w^^^ . W7^_ + («^+ . - • ("'7^+ - ■ "^'5+1 (V.9a) 

<% = {e2^(«~^)* • C(,) - e-2^(«-i)* ■ 4)} (V.9b) 

sin i9 ( * * * * ~i 




(V.9c) 



But here a nearby restriction upon the parameter b suggests itself, namely 
through the plausible demand that the Dirac currents ka flV.6[) . with their 
components being specified by ( ]V.9a|) - (1V.9c|) . must be unique (!) vector 
fields over three-space (albeit only apart from the origin r = and the z 
axis = 0,7r). Evidently this demand of uniqueness reads in terms of the 
spherical polar coordinates {r, (p} 

4(r,^?,0 + 27r) = 4(r,^9,0), (V.IO) 
and thus the values of 5 become restricted to the range 

g = lin + 1) (V.ll) 
{n = 0,±1,±2,±3,...) 

which then entails also (half- jinteger quantum numbers for the z component 
of angular momentum flV.3a)) - flV.3d[) : 

1 71 

3^ = ±(5-2) = ±2 • ^^-^^^ 

Notice here that this (half-) integrity arises as a consequence of the de- 
mand of uniqueness with respect to certain physical densities (i. e. Dirac 
current), whereas the corresponding integral quantum numbers of conven- 
tional non-relativistic quantum mechanics are mostly traced back in the text- 
books to the uniqueness requirement for the wave functions themselves (not 
the densities). The lowest values of f lV.12[) are = ±| for 5 = and 
iz = 0, ±1 for i) = ±|. Thus for the first case (S = 0) we have a fermionic 
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basis and for the second case (S = ±|) one deals with a bosonic basis. In 
this sense, a Dirac particle is said to occupy a fermionic quantum state i(j 
if the "boson number" b of its spinor basis is zero {b = 0), and a bosonic 
quantum state if the boson number b equals ±|. Observe that through this 
arrangement the fermionic or bosonic character of the quantum state of a 
Dirac particle is defined by reference to the corresponding spinor basis. 

Naturally, it must appear very tempting to interprete the results f lV.ll| - 
(]V.12|l in the sense that angular momentum Jz be quantized in units of ^. 
However, it is important to see that this "quantization" occurs in connection 
with the basis system (lV.3a[) - (]V.3d|l . not in connection with the unique wave 



amplitudes ^"'7^±, '"'iSi or the wave functions ipa, resp., which are normally 
thought to specify the angular momentum of the considered quantum state. 
Therefore one conceives here the ansatz parameter 5 to be a new quantum 
number [boson number) which does influence the magnitude of angular mo- 
mentum but does not fix it completely (i. e. the ultimate fixation, if possible, 
of angular momentum is assumed to occur through cooperation of b and the 
(unique) wave amplitudes *"'7^±(r), *'*'iS±(r)). 

2. Two-Particle Eigenvalue System 

For the two-particle systems, to be considered now, all fields which are due 
to the third particle must naturally vanish, and therefore the original eigen- 
value system for three particles becomes cut down to an eigenvalue system 
for only two particles. Thus the abstract three-particle equations flIII.40a|l - 
flIII.40c|l become converted to the corresponding eigenvalue system in terms of 
the wave amplitudes ^"'^7l±, '°1S± (a = 1,2) which are due to the ansatz flV.7a|l - 
flV.7b|) for the Pauli spinors. Each of both Dirac wave functions i/^a (a = 1, 2) 



becomes thus parametrized by four complex- valued fields ^°'^TZ±, and 
each of the two Dirac equations (]II.32a[) -( 1II.32b[) . or Pauli equations (]III.40a l 



(lIII.40b|) . resp., yields then a coupled system of four equations for these wave 
amplitudes. For the sake of brevity, it may suffice here to explicitly reproduce 
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only the case of the first particle {a — 1): 

+ - • + - • ^'^^+ + ^'^0 • (^)7^_ + [^ ^% - sin ^ • %] ■ «7e+ 

or r o<p r 

+ [i % + cos • %] ■ «5+ 1 = c • «7^_ 

(V.13a) 

^ _ 1 . ^ + h . (1)5^ + (% . + [z % + sin ^? • %] . 
or r 0(p r 



+ [i % - cos ^ • (^k^j • (^)7^+ 1 = c • 



(V.13b) 

or r 0(p r 



+ [z % + cos ^ • %] . I = Mp±Mi c . ^'>JZ, 



(V.13c) 

^ + - ■ ^ + ^ • - % . «5+ - [i % + sin ^ . %] . 

or r ocp r 

+ [i % - cos ^ . (^k^j • (^)7^_ I = + c • . 

(V.lSd) 

Since this (relativistic) system of mass eigenvalue equations for the wave 
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amplitudes looks somewhat intricate, it will be highly instructive to spend 
some effort on its detailed analysis. First let us mention that a similar sys- 
tem of eigenvalue equations does hold also for the second particle (a = 2) 
which, however, can easily be obtained by means of the particle permutation 
symmetry: 

«7^± ^ (2)7^± ^ (2)5± (V.14a) 

(^14^ = %, %} ^ (2]4M = {(%; %} (V.14b) 

Ml ^ M2 Mp ^ -Me ip^^e- (V.14c) 

Next, recall that this system (lV.13a[) -( lV.14c)) of eigenvalue equations is math- 
ematically consistent only if the boson numbers {bp, be} do obey again the 
former restriction flV.ll|l because it is only for this case that the exponential 
factors become unique functions over three-space, i. e. 

g±2j(fip-i)(/) ^ g±m</. (V.15) 

{n = 0,±1,±2,...) . 

But what is most important for the present discussion, this refers to the 
fact that the eigenvalue system flV.13a)) -f lV.14c|) can be viewed also as the set 
of variational equations due to the RST energy functional E^t] f!IV-14|) . In 



order to become convinced of this claim one merely has to represent the indi- 
vidual contributions of the energy functional in terms of the wave amplitudes 
^°'^TZ±, *°'iS± and carry through the variational procedure, which itself again 
requires the uniqueness condition ( ]V.15|) . namely in order that the neces- 
sary partial integrations can be carried out in an unambiguous way. Observe 
also that, in the most general case, all three components '"^M^, '"M^ of 
the vector potential v4a(r) will be non-zero; and this requires that in such 
a general situation the Pauli spinors ^"'(/2t(r) (lV.7a| - (lV.7b|) necessarily must 
appear as a superposition of those simple "spin-up" (~ TZ) and "spin-down" 
(~ S) components. It is only in very special situations that one can deal 
exclusively with one spin direction (either "up" or "down"), see below. 



3. One-Particle State 2p3 

2 

Surely, a single Dirac particle, moving in an external field will oc- 
cupy a fermionic state, i. e. the boson number be of a single electron must be 
adopted to be zero {be = 0). For instance, one encounters such a situation 
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for the hydrogen atom where, for the sake of simphcity, the proton mass Mp 
will be presumed to be infinite. The Dirac density ^%o(r) of such an infinitely 
heavy particle will be found to be pointlike, i. e. one puts 

= S^f) , (V.16) 

and consequently the electric potential *^%o(r) (|III.5a|) of such a point charge 
must degenerate to the Coulomb potential 

(% ^. (V.17) 

r 

The relativistic spectrum of this one-particle eigenvalue problem is exactly 
known (see e.g. ref. [IE]); and if we take this physical situation as the 
simplest test case for our RST energy functional £'[t] flIV.14[) one must expect 
to recover just that conventional spectrum [18]. Clearly, this will then entail 
some confidence into the proposed RST energy functional E[t] flIV.14)) . But 



sometimes the application of a new formalism admits also a new view upon 
the considered standard situation; and in the present case this alternative 
view provides us with the emergence of the spin degeneracy when passing 
over to the non-relativistic limit. This effect is mostly omitted in the text 
books (see, e.g., ref. |18]). 

We will treat here this phenomenon for the excited state 2p3 of a Dirac 

2 

electron in the electrostatic field (IV.17|) of a pointlike proton. Observe that 
all the magnetic fields are presumed to be zero, so that there exists locally 
no preferred direction in space for the orientation of the electronic spin. 
Therefore a "spin-up" solution must be expected to carry the same energy 
£'t as a "spin-down" solution or as any linear combination of both. However, 
this supposition is not true for the relativistic case, and the reason for this 
is the spin-orbit coupling of the angular momenta. In order to recognize this 
more clearly, one first recalls the conventional treatment of the relativistic 
hydrogen atom [18] . Here the stationary solution il'2{r) (1III.39|) of the Dirac 
equation (]II.32b|) for the electron (i. e. our second particle, a = 2) is presented 



in terms of the angular momentum eigenvectors C/'*" (lV.la| - (lV.ld[) as 

(2)(p+(f) $(r) ■ cl'^ (V.18a) 

%_{f) =^ -iW^{r) ■ cl'^ , (V.18b) 
where the radial function $(r) is given by 

$(r) = N^r'^e-'^'' . (V.19) 
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The normalization constant A^^ for this solution is 



r(3 + 2z/) 



(2/3)3+2^ 

the decay parameter (3 of the exponential function is found as 

MeC 2W 



(W = — 

as 



(V.20) 



(V.21) 



2 J 



) 



and finally the power u is given by 



■1 + 2W1 



osy 
2 J 



(V.22) 



Concerning the corresponding energy eigenvalue (-Econv, say) of this quantum 
state 2p3, it may be sufficient here to quote its non-relativistic approximation 



up to order ag [TH]: 



-etc 



1 



+ 



2-n? 211? \i + 



1 



3 

An 



(V.23) 



i. e. for the present state 2p3 with principal quantum number n = 2 and 

2 



8 128 



(V.24) 



Of course the first term of the non-relativistic expansion (]V.23[) yields the 
usual Schrodinger eigenvalue (Eg, say) of the non-relativistic hydrogen atom 



2aB n? 



(V.25) 



with the atomic energy unit (a. u.) given in terms of the elementary charge 
(e) and Bohr radius (ob) as 



Ob 



27.2116... [eV] 
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(V.26) 



But now it is especially instructive to oppose this standard treatment of 
the ideahzed hydrogen atom (1V.16[) - (1V.17[) to the corresponding RST treat- 
ment. Here it is first necessary to specify the general eigenvalue equations 
flV.14a|) -f lV.14c[) for the second particle (a = 2) which is assumed to occupy 
(as a single particle) a fermionic state, i. e. one puts the boson number §e to 
zero (§e = 0). Moreover, all the magnetic fields are also assumed to vanish 
(i. e. ^^}Ar = = = 0); and finally one introduces the real- valued wave 
amplitudes ®i?±, ^^'15'± by putting 



(V.27a) 
(V.27b) 



For such a simplified situation the system of eigenvalue equations for the 
second particle becomes cut down to the following form: 



+ mo ■ + - 

or r 



h 



(V.28a) 



or r r ov 



h 



(V.28b) 



or r 



dr 



dd 



Mo - M, 



+ - (%. 
r 



1 



r d'd 



h 



M2 - Me 

h 



(V.28c) 
(V.28d) 



This simplified system of eigenvalue equations for the electronic (real- 
valued) wave amplitudes suggests to try the following product ansatz 

^^^i?±(r, t!)) = cos (r) (V.29a) 

(2)s±(r, = sini9 ■ S±{r) , (V.29b) 

and this recasts the original eigenvalue system ( ]V.28a|) -( ]V.28d[) into the fol- 
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lowing equations for the spherically symmetric functions 

dR+(r) 



dr r h 

(V.30a) 

dS+{r) 1 n ^ ^ p ^ ^ 1 ft . X Me + M2 n , , 

ar r r n 

(V.30b) 

dR-{r) 2 u ^ ^ n ^ ^ 2 n ^ ^ M2-Me fl , , 

dr r r ft 

(V.30c) 



dr r r 

(V.SOd) 

For the solution of this set of ordinary differential equations one tries the 
nearby ansatz 

R±{r) = p± ■ r^expl-pr] (V.31a) 
S±{r) ^f±-r'' exp[-/3r] (V.Slb) 

which then yields two sets of four algebraic equations for the determination 
of the four constants p±, f± and for the ansatz parameters u and /3. The 
first set contains only the power u 

iy-p+ + 2-U + as-p- = (V.32a) 

p+ + (z/+l)-/+ + as-/- = (V.32b) 

-as ■ J5+ + (z/ + 2) ■ j9_ - 2 • /_ = (V.32c) 

-as-f+ + {iy + l)-f--p- = 0, (V.32d) 

whereas the second set relates the decay parameter (3 to the mass eigenvalue 
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p.,^^^A±Ml,.,_ (v.33a) 

h 

ff.U^'J^..}, (V.33b) 

^.,_^MlZ^,.,, (V.33C) 

,.S_^M^,.U. (V.33d) 



Naturally, the first set (IV.32a)l - (IV.32d l says that the power z/ of the RST 



ansatz (lV.31a[) - (]V.31b[) is identical to the conventional value (]V.22[) : and 



similarly the second set (]V.33a|) -( ]V.33d|) fixes the RST value of the decay 



parameter (3 also to its conventional counterpart (1V.21|) . where the mass 
eigenvalue (M2) of the electron turns out as 



with the relativistic parameter W being defined below equation flV.21[) 



This exact coincidence of the RST results and their conventional counter- 
parts mentioned above seems to signal that RST is merely a formal one-to-one 
transcription of the conventional formalism; but actually RST owns a richer 
structure. This becomes more obvious when one faces now the question of the 
"energy" carried by the considered wave function ip{r). In the conventional 
theory, this energy is quite generally identified for the presently considered 
one-particle systems with the mass eigenvalue M2C^: 

Econv = . (V.35) 

This is most easily seen in the non-relativistic limit by expanding the present 
RST result M2(? (]V.34[) with respect to the fine-structure constant as: 



M2C^ = M.c^^-^-^ ...^ (V.36) 



and comparing this to the conventional result flV.24[) . (For a general proof. 



see ref.s [T5|IT6].) However, in RST the total energy Et ( ]IV.14)I will in general 
not be found to agree with some mass eigenvalue! It is only for the present 
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oversimplified situation (1V.16[) - (1V.17[) that the RST energy can agree with 
the mass eigenvalue M2C^. This is easily seen by explicit calculation of the 
various contributions to the RST energy 

First, one makes oneself sure of the relativistic normalization condition 
(IIII.42|I for the electron by imposing upon the constants p±, f± the following 
condition: 



1 

3 ' 



r(3 + 2z/) 



(V.37a) 
(V.37b) 



Next, the Pauli spinors ^^^f±{r) flV.7a|) - flV.7b[) are inserted into the second 
mass renormalization factor 2(2) flIII.44[) which yields quite generally 

(V.38) 

i. e. in particuar for the presently considered wave amplitudes flV.29a| - 
flV.29bp and flV.31ap - flV.31bp . resp.. 



^(2) 



(2) 



2 p 

— + 2— 



(V.39) 



Here the normalization condition flV.37a|l has also been used. 

The kinetic energy 7kin{2) (HII.47P represents a more complicated problem; 
inserting here the former ansatz f ]V.7a)l -f lV.7bp for the Pauli spinors ^^±{f^ 
with the wave amplitudes ^'Tl-±^ ^S± being further specified by equations 
flV.27ap -f lV:27b|l . flV.29ap -f lV:29b]l . and flV.31ap -f lV:3Tb]l yields by straightfor- 
ward integration 



'kin(2) 



(tT^) ]^ ■ ^- " 2 ■ /- + p- ■ /+] } • (V.40) 



Obviously, both the kinetic energy Tkin(2) and the mass renormalization term 
f lV.39)l are originally built up by the integration constants of both the "pos- 
itive" /+) and the "negative" /_) type. But the relations flV.33a l- 
f lV.33dll fix their relat ive magnitudes as 

p_ = W -p^ (V.41a) 
]- = W-U- (V.41b) 
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Here the relativistic parameter W has aheady been defined below equa- 
tion (1V.21[) . but on the other hand it can be written also in terms of the 
mass eigenvalue M2 as 



W = J^^^. (V.42) 

V Me + M2 ^ ^ 

This provides us with the possibility to express both quantities by means of 
the positive-type constants alone, i. e. the mass renormalization term (]V.39|) 
and the kinetic energy (1V.40[) appear then in the following form 



Mec^ ■ Zf2) = MeC^ - '^M.c'j^ {pi + 2/^ } (V.43a) 
_ Mec^ {2WY 1 r ^ 



This form of the results will become important for the deduction of the corre- 
sponding non- relativistic limit; but for the present relativistic case one may 
use the values of the constants /+ following from the algebraic equations 
(IV.33a[) -( IV.33d[) together with the normalization condition (IV.37a[) 



/+ = -i^ ■ , ^ (V.44b) 
so that the energy contributions (lV.43a|l - (]V.43b|l adopt their final shape as 
M,c^ ■ Zj) = M,c2 - ^— ■ M,c' (V.45a) 

Thus we are left with the problem of calculating the last physical con- 
tribution to the energy £'[t] (]IV.14[) . i.e. the gauge field energy E§ f lIV3b|) 



which of course reduces to a simpler form for the two-particle systems: 

E§ = [d'fE,'E, = /dV(V(%)-(V(%) . (V.46) 
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Observe again that we need consider only the physical part E^^l^ flIV.13[) 

of the proper functional i?[T] (IIV.14|I if all the constraints are obeyed by our 
ansatz for the matter and gauge fields. But here the relativistic normalization 
constraint flIIL42|l has already been demonstrated to hold in the form flV.37a|l ; 
and furthermore it is also very instructive and satisfying to see in what 
way the electric Poisson identity (1IV.7|) is validated. To this end, one first 
computes the electric mass equivalent M^jc^ (IIII.48a)) by use of the first 
potential *^"^Mo(r) flV.17|) and the second Dirac density ^%o(0 (IV-81) 



Mgc^ = -y |d'f%„(r).%(r) 



« I H- + {"n-f + ■ (V.47) 



2 J Anr 

The straightforward calculation of this integral under use of the normalization 
constraint ( lV.37a| yields the following result 

M,§c^ = -'-^ = - ( ^) . M^c- . (V.48) 

This result is also found for the first mass equivalent, i. e. 

Mftjc^ = |dV%(f)-(%(f) = Mf^jc^ (V.49) 

which is most easily verified by simply observing the fact that the first Dirac 
density is pointhke, cf. (]V.16|l . so that the considered mass equivalent is 
essentially the value of the second potential '^Aq at the origin (r = 0) 

Mgc^ = y-%(0). (V.50) 

However, this required value of the second potential "-^o can easily be de- 
duced from the formal solution (lIII.5a[) of the corresponding Poisson equa- 
tion (IITL3]l as 

(%(0) = -asj^^Xin 

= -«ST^, (V.51) 
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cf . (1V.47[) and (1V.48[) , so that the claimed equahty (1V.49|) of both mass equiv- 
alents is immediately verified by use of equation (1V.50[) . Thus the electric 
Poisson identity flIV.7)) admits now to specify the desired gauge field energy 
E§ (roe) as 

= M(V = Mgc^ + Mgc^ = -^^M.c^ (V.52) 

But now that all three physical contributions to the total energy E'jx] 
f lIV.14|l are explicitly known, see flV.45a|) -f lV.45b[) and flV.52|) . one actually 



gets for their sum the expected numerical coincidence of the RST energy iJ-p 
and the mass eigenvalue M2C^: 

= e!^ = M,c^ ■ + 2rki„(2) + 4^ = M^c^ = E,on. ■ (V.53) 

Clearly, this is a very pleasant result in favour of our RST energy functional 
E^T] because it associates to the considered hydrogen state 2p3 exactly the 
same energy as does the conventional theory, cf. (IV. 35)1 ! On the other hand, 
the present explicit calculation of the RST energy i?T in terms of the co- 
operation and interaction of the various atomic fields yields a much more 
detailed picture (of the intra-atomic situation) than is possible in the con- 
ventional theory where the atomic energy is simply identified with the mass 
eigenvalue, without discussion of the specific interplay of matter and gauge 
fields. This feature of RST (namely to provide a more detailed picture of the 
intra-atomic situation) becomes obvious also by considering now the non- 
relativistic approximation. 



Non-Relativistic Spin Degeneracy 



The point of departure for the non-relativistic approach is again the rel- 
ativistic eigenvalue system flV.28a|l - flV.28d|l where one eliminates the wave 
amplitudes , ^^^S- in close analogy to the elimination of the Pauli spinors 
V- fllll.67a|l - fllll.67cp . i.e. one puts 



(a.), 



n I 




2MeC 1 


dr 


H I 




2MeC 1 


dr 



d 



r sin d 



+ - P>S 

r 



1 



+ 



+ 



r 



(V.54a) 
(V.54b) 
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Indeed, these approximative relations may immediately be deduced from 
their exact counterparts (lV.28a|) -( ]V.28b|) : and if this is substituted in the re- 
maining equations flV.28c[) - flV.28d|) one ends up with the desired non-relativistic 
eigenvalue equations for ^5"+ S+) and ^^^R+ R+): 



2M, 



r dr 



1 

+ ^ 



d^R^ 



+ cot^9 



dR+ 



- hSAc, ■ Ra 



2M, 



2dS. 



J.2 



sm v ov 



Es{2) ■ R+ 
(V.55a) 



-) 



(V.55b) 



Here the non-relativistic counterparts of the relativistic wave amplitudes 
^^^R+, ^^-1S'+ are denoted simply by S+; and the non-relativistic energy 
eigenvalue Es(2) is formally the same as the previous Ep(2) (1111.73b I. 

Observe also that through the elimination of the "negative" wave am- 
plitudes the "positive" wave amplitudes R+, S+ have been decou- 
pled, which is equivalent to the breaking of the spin-orbit coupling. It is 
true, both eigenvalue equations (]V.55a|) -( 1V.55b[) look very similar, but they 
are not identical because of their different angular derivative parts. But 
here one may now suppose that the angular dependence flV.29a[) -f lV.29b|) of 
the relativistic amplitudes ^^■!S'_ is left unchanged on their way to their 
non-relativistic approximations 5*+, so that the former relativistic factor- 
ization will apply also to the present non-relativistic case. And indeed, if that 
factorized ansatz (lV.29a|) - (]V.29b|) is inserted in the present non-relativistic 
equations (lV.55a[) -( ]V.55b[) one finds that both non-relativistic amplitudes 
obey the same eigenvalue equation, e. g. for the non-relativistic approxima- 
tion -R+(r) of the relativistic amplitude R+{r) 



{d^R+{r) 2di?+(r) 2 



+ -- 



2Mp I dr^ r dr 



E^2)-R+{r). 
(V.56) 



Naturally in view of the successful relativistic ansatz (lV.31a|) -( lV.31b)l . one 
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tries here for the non-relativistic situation quite analogously 



R+{r) = po ■ r"" exp[-(3or] (V.57a) 
5+(r) = /o ■ exp[-/9or] . (V.57b) 

The ansatz parameters and (3^ are then determined from equations (IV.56[) 



as 



Uo = 1 (V.58b) 

which of course is nothing else than the non-relativistic limit of their orig- 
inal counterparts (3 ( ]V.21|) and u (1V.22[) . Furthermore, the non-relativistic 



eigenvalue £s(2) is found as 



£k(2) = 'l4-Mec' = (V.59) 



which again is the non-relativistic limit of the mass eigenvalue M2(? f lV.34[) 



or the convential energy £tonv ( 1V.24|) . (Recall here that the non-relativistic 
limits are deduced from their relativistic counterparts through expanding the 
latter with respect to the fine-structure constant as-) 

But concerning the non-relativistic form of the wave function itself, one 
is referred back to the original ansatz flV.7a|l for the "positive" Pauli spinor 



^\)+{f) and observes that the corresponding "negative" spinor ^^\p_{r) flV.7b[) 



is neglected for the non-relativistic approach. Thus the non-relativistic spinor 
solution for the present fermionic (5e = 0) state 2p3 reads 



(2)^^(f) cos^ R+{r) ■ Q'^ +e''^ sin dS+{r) ■ Q' ^ 

= re-^°"|j9oCos^?Co^'^ +/oe^'^sint?Co^'~^| . (V.60) 

Imposing now upon this solution the non-relativistic normalization condition 
analogous to (]III.100|) yields the following constraint for the non-relativistic 
normalization constants Po, fo'- 

pI + 21' = -L. (V.61) 
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In order to satisfy this condition trivially, one introduces the degeneracy 
parameter B through putting 

Po = cos e (V.62a) 

/o = -^sine, (V.62b) 

which then recasts the non-relativistic solution (1V.60| into its final form 

%+{r) ^-^j^lcosQcos'd-Cp + . (V.63) 

V8a| I V2 J 

This non-relativistic RST result should now be compared to the non- 
relativistic limit of the conventional form flV.18a|) of the relativistic 2p3 state. 
To this end, one first expands the relativistic normalization factor A^^ (IV.20|) 
with respect to the fine-structure constant as which yields in the lowest order 
of approximation 



and if furthermore the decomposition of the basis spinor d"^ with respect 

11 1 

to our fermionic (5e = 0) basis {Co*'^ , Co^' ^ } is applied, i. e. 

= V2|cos^?-Co^'^ -^e^'^sin^9-Co^'"^| , (V.65) 

then the non-relativistic limit of the conventional Pauli spinor (1V.18a)l ap- 
pears as 

r 

^ j COS 7? ■ Co^'^ - - e''^ sin d ■ Q''^ \ . (V.66) 

Thus the comparison of this conventional result to the degenerate RST 
state (]V.63|) fixes the degeneracy parameter 9 to 

cos e = yi (V.67a) 
sinG = -y|. (V.67b) 
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So we see that the conventional approach provides a somewhat too sim- 
ple transition (1V.66|) to the non-relativistic limit by unnecessarily fixing the 
degeneracy parameter which, however, by the RST approach flV.63|l is 
correctly left undetermined! 

Finally, we have to reassure that the non-relativistic RST functional 
E[T] flIV.21|) actually does confirm the degeneracy effect, i. e. its value upon 
the non-relativistic RST solution flV.63[) must be independent of the degener- 
acy parameter 9 introduced by equations (]V.62a|) -( 1V.62b[) . In order to verify 
this by a few brief arguments, one first observes that through the present ne- 
glection of the exchange effects and of the magnetism that non-relativistic 
functional E[t] (IIV.21)I becomes cut down to the following form 

E[T] e[?] = ^kin + 4^ + Ap(2)-ND(2) + Ag-N|. (V.68) 

But here the non-relativistic normalization condition flIII.131|l is obviously 
satisfied by our present result ( ]V.63|) . and the same is true also for the non- 
relativistic Poisson identity ( ]IV.20a)l . This is easily seen by explicitly cal- 
culating both mass equivalents Mj^^jC^ and M^jC^ (]III.48a| which yields by 
simple integration 

Mg = Mg = , (V.69) 
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so that the gauge field energy becomes 



E§ = M^c' + m|c^ = -— . (V.70) 

Observe here that even this gauge field energy alone is independent of the 
degeneracy parameter 0. But this then must hold also for the kinetic 
energy -Ekm(2) of the second particle. Indeed, the explicit integration for 
-£i;in(2) fllll.llTal) with reference to the present solution (^l/j+(r) flV.63|l im- 
mediately yields 

iii„(2) = (V.71) 
Thus adding up both contributions (]V.70|) and (1V.71[) to the total non- 



relativistic energy e[^, frTOI) yields 



< = = -^«i^ec^ (V.72) 
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in perfect agreement with the short-hand expansion of the relativistic result 
(1V.35| - (1V.36[) . Observe also that in the non- relativistic limit the total energy 
e[^j f l\r72|) coincides again with the non-relativistic eigenvalue -Es(2) f !V-59|) ! 

Summarizing, the constructed RST energy functional for the fermionic 
one-particle systems meets with all the logical demands which one can rea- 
sonably impose on such a functional. Surely, this yields sufficient motivation 
to test now also its viability for the more complicated physical situations, 
such as, e.g., a non-trivial two-particle system where the Dirac particles 
occupy bosonic states! 



4' Simplest Bosonic State (5 = |) 

A bosonic state naturally occurs for the positronium groundstate [TT]. 
For such a bound system of a positron (a = 1) and an electron (a = 2), 
the former two-particle eigenvalue system flV.13a|l -f lV.14c|) does also apply, 
however with identical rest masses and mass eigenvalues, i. e. 



M, = Mp = M (V.73a) 

Ml = -M2 = -M, (V.73b) 
1 
2 



be=bp = - . (V.73c) 



Furthermore, it is well known that positronium emits a one-particle spec- 
trum, i. e. the two oppositely charged constituents must occupy (almost) 
identical quantum states. This requires that the Dirac charge densities ^"'^ko{r) 
are identical 

%(f) = (%(f) = («A;o(f) (V.74) 

so that the Maxwell densities ^"'%{f) fill. 43a|l - fill. 43b|l differ in sign: 

%{f) ^ ~%{r) = ^%{r) , (V.75) 

which then also holds for the electric potentials *'*]4o(r) generated by these 
densities according to the Poisson equations fllll.3|) : 

%(f) = -%(f) = »Mo(f) . (V.76) 

Concerning now the magnetic fields Ha{f), one encounters two possibilities: 
either the magnetic fields are antiparallel {Hi{r) = —H2{r), para-positronium) 
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or they are parallel {ortho-positronium [11]): 



Ai{r) = A2{f) 
HAr) = mf] 



Mr) 



(V.77a) 
(V.77b) 



(In the conventional theory, these are the para-states """Sq.) Clearly, such a 
coincidence of the electromagnetic fields requires a similar coincidence of the 
corresponding Dirac three-currents kaif) (or Maxwell currents jai'f^), resp.). 
Indeed, it is easy to see by reference to the Poisson equations (1III.4| that 
the magnetic demands (lV.77a| -f lV.77b|) imply the identity of the Maxwell 
currents jair) or antiparallelity of the Dirac currents ka{r), resp. 



3i[r) =j2[r, 



ki{r) = -k2{r) = kp{f) 



(V.78a) 
(V.78b) 



Of course, such a far-reaching coincidence of the matter densities and 
gauge fields of both particles must be induced by the corresponding rela- 
tionship of the wave functions. But here it is a simple thing to satisfy both 
density requirements flV.74[) and flV.78b[) : indeed, resorting to the original def- 
initions of those densities in terms of the Pauli spinors ^°'^if±{f), cf. ( IIII.123|) 
and 0IV.22|) . one puts 



(V(r) 



(1), 



(V.79a) 
(V.79b) 



which then reads in terms of the wave amplitudes (]V.7a[) - (]V.7b l 



«7^+ 


= ®7^+ 


= *?R+ 


(V.SOa) 




= = 




(V.SOb) 








(V.SOc) 




= 


= '^IS 


(V.SOd) 



Accordingly, the relativistic two-particle eigenvalue equations (lIII.40a[) - (lIII.40b[) 
must collapse to the following one-particle eigenvalue problem: 



h 



(V.81) 
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This abstract system may be rewritten also in terms of the (real-valued) 
wave amplitudes *??±, ^^^S±, namely by simply introducing those identifi- 
cations flV.80a|) -f lV.80d|) into the two-particle system flV.13a)) -f lV.14c|) which 
then collapses to the following one-particle system 



dR+ 1 dS+ 
+ 

dr r dd 


- *]4o ■ R- = 


M + M, 


R_ 


(V.82a) 


dS+ 1 dR+ 
dr r dd 


- *]4o ■ S. = 


M + M, 


S^ 


(V.82b) 


ld{r-R_) Ids. 
r dr r d-d 


+ ■ R+ = 




R+ 


(V.82c) 


ld{r-S_) IdR. 
r dr r dd 


+ *Mo ■ S+ = 






(V.82d) 


magnetic field is ne 


glected {Ap =^ 


0), see ref. 


m 


for its inclusion, 



and the modified wave amplitudes R±, S± are defined through 



R± = Vrsind ■ *!R± (V.83a) 
S± = Vrsin^ ■ »!S± . (V.83b) 

Naturally, this one-particle eigenvalue system (]V.82a|) -( ]V.82d[) couples to 
merely a single gauge potential *Mo, cf. (1V.76|) : and consequently the two 



Poisson equations (1III.3| become contracted to a single one for this residual 
potential 

62 _|_ c2 _|_ 52 I c2 

A*Mo = -47ras'% = -as ^+ + + ^- + (v.84) 

r smv 

whose formal standard solution is as usual 

'%(r^) - /■ Rl{r',^') + Sl{r',^') + Rl{r',^') + Sl{r',^') 



4:71.1 r'sint?' ||r — r'|| 

(V.85) 



Comparing the present bosonic system flV.82a|l - flV.82d|) together with the 



Poisson equation flV.84|) to its fermionic counterpart flV.28a|) -f lV.28d|l . the 



crucial difference must refer of course to the electric potential: whereas the 
Coulomb potential '■^o (1V.17|) is fixed from the outside for the fermionic case 
(so that one could obtain an exact solution), the present interaction potential 
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*]4o ( 1V.85|) is part of the dynamics and therefore must be determined simul- 
taneously with the wave amplitudes. Of course, this circumstance spoils the 
possibility of finding an exact solution for the present bosonic case; and thus 
one has to resort to some approximative method, e. g. the variational method 
based upon just the energy functional E^t]- 

In this sense, it becomes necessary to specialize also the general result 
for E[T] ( ]IV.14)l to the present situation where the magnetic and exchange 
fields are neglected. Observing here first the circumstance that the identifica- 
tions (lV.80a[) -( lV.80d[) ultimately leave us with the one-particle fields R±, S± 
f lV.83a[) -( ]V.83b[) . the energy functional E^t] is found to be simplified to the 
corresponding one-particle contributions, i. e. 

^ 2Mc2 ■ J2 + 4f^.^ + E§ + (Ad(i) + Ad(2)) ■ Nd + A? ■ ^^^"^ 



El 



[T] 



(V.86) 



with the one-particle objects being given in a self-evident way through (d^r == 
r drd'd) 



7-2 _ 72 _ 

^(1) - ^(2) - 
Tkin(l) = 7kin(2) 



(V.87a) 



kin 



A2^ J 6 dR+ R+d{rR_) ~ dS+ S+d{rS^ 



Id'r <! R_ 



1 

+ - 

r 



dr 



r dr 



~ dR^ 



dr 




dS+' 






1 


(V.87b) 



'D(l) 
(e) 



Nr 



E. 



R 



(e) 



E, 



'D(2) -. 

he 

47rac 



dV- \Rl + Si + R^ + S^_ 



d¥||V*]4nll^ 



Rl + Sl + Rl + 



. 



(V.87c) 

(V.87d) 

(V.87e) 
(V.87f) 



But with these arrangements it becomes an easy exercise to deduce both 
the one-particle eigenvalue equations flV.82a|l -f lV.82d|) and the Poisson equa- 
tion flV.84| as the variational equations due to the truncated RST energy 
functional ^[t] flV:86|) . 
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5. Positronium Groundstate (I^Sq) 



The bosonic variational principle flV.86[) with the corresponding varia- 



tional equations flV.82a[) - flV.85[) describes that part of the positronium spec- 
trum which is due to the boson number 5e = 5p = | and which additionally 
admits the use of real- valued wave amplitudes *R±, ^^^S± f lV.80a| -f lV.80d|l . 
The question of completeness of the positronium spectrum will require an 
extra discussion, but the groundstate (as the state of highest symmetry) 
will surely be a member of the present subset defined by the energy func- 
tional (]V.86[) . As a test of this supposition one may assume that the ground- 
state wave amplitudes R±{r,'{}), S±{r,'d) are approximately spherically sym- 
metric, i. e. we put for the groundstate 

^^(r,i?) ^ i?±(r) (V.88a) 
S±{r,d) ^ S±{r) , (V.88b) 

and simultaneously we pass over to the non-relativistic approximation. If 
RST has something to do with reality, one must recover (at least approx- 
imately) for this situation the non-relativistic positronium groundstate to- 
gether with the excited states, because of the non-relativistic angular momen- 
tum degeneracy. Now in order to verify this expectation, one first looks for 
the non-relativistic version (Ef°lj) of the original energy functional i?[T] f lV.86|) 
in the spherically symmmetric approximation flV.88a[) - flV.88b[) . i.e. 



E;?] = 2^kin + ^^f + 2AsNd + A§n|, (V.89) 

with the non-relativistic version E'kin of the one-particle kinetic energy Tkin 
( IV.87b[) being given by 




Tkin ^ i?kin = TTTj-- dr r — ±^ , (V.90) 



furthermore the spherically symmetric form of the gauge field energy E^ 
flV.87d)) is (*]4o(r,^?) Ao{r)) 




(V.91) 
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the mass equivalent M^^c^ (lV.87e)l becomes 

oo 

M(V ^ M^V = -^c-^ jdrrAoir)Rlir) , (V.92) 







and finally the normalization condition flV.87c|) appears in the spherically 
symmetric approximations as 

oo 

Nd ^ Nd = ^ jdrrRlir) -1 = 0. (V.93) 







The non-relativistic form of the eigenvalue equations flV.82a|) - flV.82d[) 



reads now for the spherically symmetric wave amplitude R+{r) 
^2 1 d / ^ _ dR+(r) _ ^^^^^^^ ^ _ ^^^^^ 



2M r dr \ dr 

Of course, this equation can also be deduced by extremalizing the energy 
functional Ej^j (]V.89|l with respect to the wave amplitude R+{r) (5e[Jj = 0). 
For the sake of simplicity, the non-relativistic spin degeneracy is suppressed 
here (-w <S'+(r) 0). Finally, the Poisson equation for the spherically sym- 
metric gauge potential AqIt) appears as the corresponding variational equa- 



tion of Eg as 



Id, dAo{r)\ ^ . TT Rlirj ^^^^ , 



The present positronium groundstate problem (]V.89|) - (]V.95|l has been 
treated in great detail in a preceding paper, see equations (III.27)-(III.34) 
of ref. [lOj; and it was found that the corresponding RST predictions for the 
non-relativistic positronium spectrum are in acceptable agreement with the 
conventional spectrum 

p2 1 fi 8029 

K = 0,l,2,...) 
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see fig. 2 of ref. [TU]. It is true, there are certain deviations (~ 10%) of the 
RST predictions from the conventional results (1V.96[) but these deviations 
can be attributed to just the use of the spherically symmetric approxima- 
tion. This logical situation lends itself now to a further test of the RST 
energy functional flV.86|l : namely one can apply the (non-relativistic) varia- 
tional principle 5'E[t] = (due to the non-relativistic approximation ']E[t] 
of Et ( IV. 86)1 ) to that excited positronium state (i.e. 2-'-Pi) which has in 
the spherically symmetric approximation the same orbital symmetry as the 
fermionic one-particle state 2p3 treated in the preceding subsection. Because 

2 

of the non-relativistic angular momentum degeneracy, the non-relativistic 
RST energy 'Et of that excited positronium state 2^Pi (conventional classi- 
fication) should turn out according to ( 1V.96|) as {up = 1) 



'Ex 



e2 1 e2 



2iPi 4aB (1 + 1)2 16aB 



-1.70... [eV] . (V.97) 



Of course, we are not able here to verify this claim flV.97| exactly, because 
the coupled RST system of eigenvalue and Poisson equations is not exactly 
solvable, not even in the non-relativistic limit; but one could accept the 
verification of (]V.97[) as a support of the RST energy functional if the cor- 
responding RST prediction turns out to meet with the requirement (1V.97|) 
up to that magnitude of inaccuracy (~ 10%) which showed up also for the 
groundstate (see fig. 2 of ref. [lOj) and is to be attributed to the use of the 
spherically symmetric approximation. 

6. Excited Positronium State (2^Pi) 

Concerning the excited positronium states, it must appear as a matter 
of course that exact (relativistic or non-relativistic) solutions of the coupled 
matter and gauge field system (]V.82a|) -( ]V.85[) will not be attainable. There- 
fore we have to resort again to some approximation which we adopt here in 
a two-fold way, namely in form of a combination of the non-relativistic limit 
and the spherically symmetric approximation. 

Turning first to the latter kind of approximation, one can of course not 
take over the former assumptions (IV.88a)l -f lV.88b|) for the groundstate to the 
presently considered excited state; but rather for the wave amplitudes R±, S± 



we resort to that angular dependence flV.29a)) - flV.29b[) of the fermionic one- 
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particle state 2p3, i.e. we put 

(r, 'R± (r) ■ cos ( V.98a) 

S±{r,^) 'S±{r) ■ sini? . (V.98b) 

This assumption for the angular dependence of the wave amplitudes is further 
complemented by the assumption of strict S0{3) symmetry for the gauge 
potential 'Pi4o(r, ■!?), i.e. we put 

(mo(r,t?) Ao(r) . (V.99) 

These symmetry assumptions recast the exact positronium eigenvalue prob- 
lem (IV.82a|l -f lV.84)l to the following approximate (but still relativistic) form: 



^^f^ + --S+(r)-Ao(r).'i?_(r 
dr 



dr r 

1 A (r'R_(r)) - 1 . S_(r) + Ao(r) • 'R+{r 



r 



r 



M + M, 

h 

M + M, 

h 

M-M, 

h 

M-M. 



c-'R^{r) (V.lOOa) 



c-'S^{r) (V.lOOb) 



c-'R+{r) (V.lOOc) 



h 



--as 



c-'S+{r) (V.lOOd) 



nl + sl + 'R^_ + si 



-^{r'S4r))--.'R4r) + Ao{r).'S+{ 

(V.lOOe) 

Here it is also interesting to point at the differences of the present two-particle 
eigenvalue system ( ]V.100a|) -( lV.100d|l compared to the former one-particle 
eigenvalue problem (1V.30a|) - (1V.30d[) . Clearly these differences concern the 
interaction potential Ao(r) and are responsible for the different energy spectra 
of the one- and two-particle cases, resp. Recall here the fact that in the non- 
relativistic conventional theory these differences simply consist in replacing 
the rest mass M of the one-particle problem (hydrogen) by ^ order to 
obtain the two-particle spectrum flV.96| ! 
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It is also important to observe that the assumption of angular dependence 
(lV.98a[) -( ]V.98b[) together with the strict spherical symmetry of the gauge po- 
tential Ao(r) flV.99[) does not spoil the principle of minimal energy. Indeed 
it is an easy exercise to demonstrate that the present approximate eigen- 
value equations flV.100a|) - flV.100e[) are the variational equations due to the 
functional 'E[t] which arises from the original functional i?[T] f!V-86[) through 
applying the assumptions (IV.98a| -f lV.99)) of spherical symmetry: 



El 



[T] 



'^[T] = 2Mc2.'zV4'rki„ + 'i^g+(AD(i) + AD(2))-% + AS^-'iV^^ 



(V.lOl) 



Here the relativistic mass renormalization 'Z'^ is deduced from the original 
frvrSTa)) as 



'Z' 



— I drr 



{'R+ir)Y + (S+(r))^ - {'R.ir)Y - CS-ir))'] , (V.102) 



next the original kinetic energy Ikin flV.87b|) transcribes to 'T^in as 

oo 

d'R+{r) 'R+{r) d{rR_{r)) 



ikin 



TT 



kin 



— he drr< R_ (r) ■ 



dr 



dr 



+'S(r) 



dS+{r) S+{r) d(r^_(r)) 2 



dr 



dr 



+ -[R+{r)-'S4r)+'S+{r)-n^{r)] 



(V.103) 



furthermore the gauge field energy f lV.87d[) becomes a purely radial in 
tegral due to the symmetry assumption flV.99 l 



E, 



(e) 
R 



'E 



oo 

/ dr r^ 

tts J 





dAo(r) 
dr 



(V.104) 



the normalization constraint Nd ( ]V.87c|) appears now as 

oo 

Nb ^ 'Nn = ^y"drr{('i?+(r))'+(S+(r))'+('i?_(r))'+(»_(r))'} 



= 0, 
(V.105) 
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and finally the electric Poisson identity (lV.87f|l reads 



iV|) ^ 'N§ = 'E§-'M%'' (V.106) 
with the electric mass equivalent M^'^^c^ (lV.87e[) adopting its new form as 



M(V ^ 'M^'^c^ = --he |drrAo{r)\{'R+{r)Y+(;S+{r)Y+{'R4r))'' 



+ CS-ir))' 

(V.107) 

The next step of approximation concerns again the neglection of the rel- 
ativistic effects. Here, the procedure is as usual: Calculating approximately 
the negative Pauli amplitudes (r) and S*- (r) from the equations (]V.100a|l - 
(IV.innbjl as 



and substituting this into the residual equations (]V.100c|) -( ]V.100d[) yields 
again the non-relativistic eigenvalue equations for the positive Pauli ampli- 
tudes 5"+, e.g. for -R+(r) as the non-relativistic version of R+{r) 

+ -^cMr)-R+{r) = E,-R+{r) 



2M I dr^ r dr 




(V.109) 



The same equation is also found for the spin-down component S+{r), but 
one may renounce on this if one is willing to omit the corresponding spin 
degeneracy. The non-relativistic eigenvalue equation flV.lOQjl must be com- 
plemented by the non-relativistic version of the Poisson equation flV.84[) . i.e. 



by neglection of the negative Pauli amplitudes 'S- and by omission of 
the positive spin-down amplitude 5*+ (due to the spin degeneracy): 




,2 



r 
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Clearly, the coupled non-relativistic system (1V.109[) -( 1V.110[) can again be 
understood to represent the variational equations due to an associated en- 
ergy functional (Ej^j, say) which deduces from its relativistic predecessor 

'E'[T] (1V.101|) just through the approximations (IV. 108a[) - (1 V. lOSbf . namely 

%] ^ e(?, = 2^kin + '^|^ + 2As-ND + Ag-N| (V.Ul) 



with 



E..-^U^rrl(m]\^, ,V.U2a) 



2M4 7 1 I dr / r2 



oo 



Nd = - jdrrRX{r) -1 = (V.112b) 







n| = 'Ef - m(^)c2 = (V.112c) 

oo 

M(V = -hc^ jdrrAo{r){R+{r)Y . (V.112d) 



The Lagrangean multiplier As is half the non-relativistic limit of the sum 

Ad(i) + Ad(2), namely 

Ad(i) + Ad(2) = (Ml - M2)c^ i-Mc^ - E,) - (Mc^ + E,) -2E, = 2As , 

(V.113) 



cf. flIII.51ap - flIII.51bp . 



The non- relativistic form eJ^j f lV.llip of our RST principle of minimal 
energy is now just the right point of departure for calculating approximately 
the (non-relativistic) energy of the considered positronium state 2^Pi. How- 
ever, concerning the intra-atomic situation from the viewpoint of various 
energy contributions, this variational principle ^Ej^^j owns some interesting 
features to be first elaborated in some detail. 



Virial Theorem 

The total energy Ej^ flV.lll|) of any solution R+{r), Ao(r) of the eigen- 
value problem flV.109[) - flV.110| is obviously built up exclusively by the sum of 
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the kinetic energy E'kin and the gauge field energy 'E^, provided the consid- 
ered solutions satisfy the normalization constraint flV.112b[) and the Poisson 
identity flV.112c| . Thus one first has two equations for the four energy con- 
tributions i?*, E'kin, '-Er, M^'^^c^; namely 

E^^ = 2^kin + '-Ef (V.114a) 
'e|^=M(^)c2. (V.114b) 

Additionally, a third relation can be obtained from the eigenvalue equa- 
tion fIV. 109)1 by multiplying through that equation with -R+(r) and integrat- 
ing over, using the normalization constraint flV.112b[) : 

= Ekin + M^V. (V.115) 

But the interesting point here is now that a fourth relation does exist, i.e. 

'E§ = 4- (M^V-S,) , (V.116) 

and this fact then admits to express all energy contributions in terms of the 
electric mass equivalent M^^^c^ 

= ^M(V (V.llTa) 

E^) = iM(V = ^E, (V.llTb) 

^kin = -^M(V. (V.117c) 

The total kinetic energy of the two-particle system is thus revealed as just 
half of its electric interaction energy (up to sign), i. e. 

2^kin = -^M(V = (V.118) 

and furthermore the total energy equals minus the total kinetic energy 
(i. e. twice the kinetic energy per particle) 

E^^ = -2^km = ^M^V. (V.119) 

This result is the well-known virial theorem of classical and quantum me- 
chanics [T9| |20] which in those conventional theories does apply to systems 
of particles being bound together via the Coulomb force. 
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In order to verify the vahdity of this virial theorem for the present RST 
situation, one merely has to supply the proof of the fourth energy rela- 
tion flV.116|) . For this purpose, one reconsiders the present variational prin- 
ciple ^Ej^j = (]V.111|) and subjects the corresponding extremal field config- 
urations to a global spatial dilatation (r ^ A • r, A = const), i.e. 



io(r) 



i?+(Ar) 
io(Ar) 



(V.120a) 
(V.120b) 



According to this common modification of the matter and gauge fields, the 
various contributions to the energy functional Ej^j (IV.lllf undergo the fol- 
lowing changes: 



Nd 



E'kin , no change 
A ^ 

^ (Nd + 1) - 1 

^ ■ m(V . 
A"' 



(V.121a) 
(V.121b) 

(V.121c) 

(V.121d) 



Putting these individual changes together, the energy functional E^'' reacts 
as follows to the dilatational transformation (IV. 120a|) - (] V. 1 20b[) 



E, 



E^^(A) 



2-Eki 



kin 



^ (Nd + 1) - 1 



+ Ag 



(V.122) 



But since we started from an extremal field configuration (5E|^j = 0), we 
obviously must have 



dE?^(A) 



dA 



0, 



(V.123) 



A=l 



and this just reproduces the desired fourth relation flV.116|l . since the La- 
grangean multiplier Ag is to be identified with the non-relativistic energy 
eigenvalue E^, see equation flV.113|) . 
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In the present context of this RST virial theorem, it is important to re- 
mark that the extremal field configuration must not necessarily be an exact 
solution of the associated eigenvalue problem flV.109)) -f lV.110|) . Rather, it is 
merely necessary that the dilatated fields flV.120a|l - flV.120b[) remain members 



of that subspace of trial functions which is considered for extremalizing the 
energy functional E^^ flV.llip ! This then means that the RST virial theo- 
rem does apply also to the approximate variational solutions of the original 
eigenvalue problem (1V.109|) - (]V.110[) . In order to present an example of this 
type, consider the set of hydrogen-like wave functions which appear as the 
product of some polynomial Pu{r) times an exponential function 

R+{r) ^ N ■P^{r)- exp[-/?r] (V.124a) 
P,{r) = ^anr'' (V.124b) 

n=0 

(A^ = normalization constant) . 

Since the polynomial coefficients and the decay parameter (3 are un- 
derstood to work as the variational parameters for extremalizing the func- 
tional E^J'' flV.lll|l . the dilatational transformations flV.120a|) obviously do 



not change the hydrogen-like character of the trial functions, since those 
transformations merely induce some irrelevant change of the variational pa- 
rameters {a„,/3}. The virial theorem (1V.118|I - (]V.119[) must then also hold 



for the extremal configurations within the subset of hydrogen-like wave func- 
tions, see the treatment of the groundstate in ref. [10]. 



Energy of Excited State 2^Pi 



Those trial functions of the type (IV.124ap - (IV.124b l are well-suited in 



order to estimate the (non-relativistic) excitation energy of the considered 
singlet state 2^Pi (this is the conventional classification; the corresponding 
RST classification requires an extra discussion). However, in order to keep 
the calculations as simple as possible, we resort to the simplest possible form 
of trial functions by regarding only one term of the polynomial Py flV.124b|l . 
i. e. we try 

R^{r) = Nr'^e-^''. (V.125) 

It is true, this is the same form as for the the precedent one-particle situation 
(|V.57a[) - flV.57b[) : but of course the variational parameters u and (] must adopt 
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here other values than those being specified by equations (lV.58a| -f lV.58b l. 



Despite such a simple form of the trial functions, the subsequent calculations 
represent a brief demonstration of how to make use of the RST principle of 
minimal energy. 

In the first step, one calculates the normalization constant in terms of 
the variational parameters u and (3, which yields by reference to the normal- 
ization constraint flV.112b l 



= i.S^^)^, (V.126) 

71 r(2// + 2) ^ ' 

Using always this value for the normalization constant A^, one can omit the 
constraint flV.112b[) in the energy functional flV.lllf . And once a normalized 



trial function R+{r) is at hand now, one can immediately compute the kinetic 
energy E]^[^ (]V.112a)l in terms of the variational parameters. The general form 



looks as follows 

^kin = ^m'-eUv) (V.127) 

where the dimensionless function £kin(i^) of the variational parameter v is 
given for that simple trial function fIV. 125)1 by 

provided is a (half-)integer: v = 0, |, 1, |, 2, |, ... (for the case of more 
complicated functions £kin see equation (IV.22) of ref. [lOj). 

In the next step we have to determine the gauge potential AQ{r) from the 
Poisson equation (1V.110[) . because only in this case we can omit the Poisson 
constraint flV.112c[) in the energy functional Ej^j flV.lll[) . Strictly speaking, 
it is not even necessary to determine the gauge potential Ao{r) itself but 
rather it suffices to determine its field strength since only this enters 

the physical part of the energy functional e|^j (provided the Poisson identity 
is satisfied). The general form of the field strength is chosen to look as follows 

m 



dr 



|l_e-2^^[l + /(r)]} , (V.129) 
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and if this is substituted into the Poisson equation (1V.110[) one gets as the 
differential equation for the ansatz function /(r) 

Mr) _ 2/3 [1 + /(r)] = --iVV^+^ 



dr L •> y n 4 

(2/?r)2^+i . (V.130) 



r(2z/ + 2) 

The solution of such a simple differential equation can easily be deter- 
mined by first introducing the dimensionless variable y through 

y = 2(3r (V.131) 

which recasts the original differential equation into the following form 

^lM-\l + f(y)] = (V.132) 

dy L ■'^^'^ r(2/y + 2) ^ ' 

with the corresponding solution 



n=l 



provided the variational parameter z/ is a (half-)integer: z/ = 0, |, 1, |, ... Since 
furthermore the field strength ( IV. 129)1 reads in terms of the dimensionless 
variable y 



2 



, , i^PasY^ ^, (V.134) 

dy J ?/4 

the gauge field energy 'E^ ( 1V.104[) adopts the following form 

'^1^ = -as^c-2/5-epot(//) . (V.135) 

Here the dimensionless function Spoti^) of the variational parameter u is given 
by 



-pot I 



j,, ji^fZ^y , (V,36) 
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i. e. for our special case (IV. 133)1 

, - _ (n + m-2)! 1 



n! m! 2"+"'-i ^ n(n - 1) 



n,m=l n=2 



(V.137) 

Observe here that both functions ekinl^^) f!V-128|l and £pot('^) f!V-137|) are de- 



fined here exclusively for (half-) integer values of the variational parameter 
V = 0,|,1,|,..., since for general ly such integrals as showing up in equa- 
tion flV.136| could not be given in a simple analytic form! But on principle, 
1^ is to be considered as a continuous variable for which both functions £km('^) 
and £:pot(i^) do surely exist. 

But now that all constraints are satisfied and both the kinetic energy 
-^kin (1V.127| and the gauge field energy 'E^ (1V.135[) are explicitly known in 
terms of the variational parameters u and (3, one can build up by means of 
them the total energy Ei^^ (IV.lll)! 

F.^^\(3,u) = 2Eki„ + '^|^ = e^aB (2/3)2. ekinW- 6^(2/3) -SpotM (V.138) 
(ob == ^ ^ ■ ■ ■ Bohr radius) . 

According to the principle of minimal energy, the minimal value of this 
function ]e|J^(/3, p) yields the energy of the considered state in both the 
non-relativistic limit and in the spherically symmetric approximation. Turn- 
ing first to the minimalization process with respect to the decay parameter 
/9, one has 

= ^^!^ = 2e2[4aB-6ki„(i^)-/3-epot(^)] , (V.139) 
which fixes the minimalizing value 13^ of the decay constant jS to 

^ = T^-^- (V.140) 

Now substituting this back into the kinetic and gauge field energies i?kin f!V.127p 
and 'E^^ (jVHSj one arrives at 

{2(3.r ■ sU^) = ^ ■ (V.Mla) 



E, 



kin 



P* 2M SOb £km(i^) 

2 c-2 



-e^ {2(3,) ■ e^,,{u) = -— ■ , (V.141b) 
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and this result just realizes the two-particle virial theorem (1V.118| ! (Clearly, 
the one-particle virial theorem is also satisfied, cf. equations (1V.70| - (1V.72)) ). 
Accordingly, the total energy becomes 

= P£p^ ^ . s{u) (V.142) 

16aB £kin(l^) 



where -Etonv is the prediction (1V.96[) of standard quantum mechanics for prin- 



cipal quantum number rip = 1; and the function S{i') 

S{u) = (V.143) 

contributes to Ej!^ by its maximally possible value (according to the prin- 
ciple of minimal energy). 

If this maximal value of the function S{iy) turned out to be unity 
(i.e. S"* = 1), the RST prediction E^'' f lV.142|) would exactly agree with 
the conventional prediction i?conv|„ _^ (]V.97|l . Perhaps this hypothesis of the 
"exact" coincidence of the RST and conventional predictions is correct for 
the true non-relativistic limit of RST, which does not rely on the spherically 
symmetric approximation flV.98a|) -f lV.99[) of the gauge potential ^P]4o(r, t?). 
But as the discussion of the groundstate l"^So in ref.s [TOl [11] has shown, the 
latter type of approximation induces a deviation of magnitude ~ 10%, which 
is then to be expected also for the present case of 2^Pi. Indeed, approximating 
the properly continuous function S^iy) (1V.143[) (with the discrete values as 
shown in the figure below) by an interpolating polynomial of 11th order 
yields a maximum = 0.9119... at z/* = 1.7942... and thus realizes the 10% 
expectations (see fig.). 
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Energy Level Diagram 
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Minimal value of the energy E^p^ ({V.14^ 

The coincidence of the RST prediction flV.142[) with the conventional 
prediction £i;onv|„ f lV.97[) would occur if the maximal value of the func- 
tion S{iy) (1V.143[) were unity {S{iy) =^ 1). Since, however, the maximal 
value 5* due to the trial function (]V.125[) is found as 5"* = 0.9119..., one 
ends up with a deviation of (roughly) 10%, i. e. the RST prediction for 
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the energy is —1.55... [eV] as compared to the conventional prediction 
£^conv|„ _^ = 1-70... [eV]. This deviation is to be ascribed partly to the use of 
the spherically symmetric approximation and partly to the use of an oversim- 
plified trial function (1V.125[) . Thus it remains to be clarified to what extent 
the RST and conventional predictions would agree if better approximation 
techniques for the RST calculations could be applied. 
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